Carving Out the Space of 4D CFTs 



David Poland^'^, David Simmons-Duffin^, and Alessandro Vichi^'^ 

^Jefferson Physical Laboratory, Harvard University, Cambridge, MA 02138, USA 
"^School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA 
^ Institut de Theorie des Phenomenes Physiques, EPFL, Lausanne, Switzerland 
^ Ernest Orlando Lawrence Berkeley National Laboratory, 
University of California, Berkeley, CA 94720, USA 

September 27, 2011 

Abstract 

We introduce a new numerical algorithm based on semidefinite programming to effi- 
ciently compute bounds on operator dimensions, central charges, and OPE coefficients 
in 4D conformal and M = 1 superconformal field theories. Using our algorithm, we 
dramatically improve previous bounds on a number of CFT quantities, particularly for 
theories with global symmetries. In the case of S0(4) or SU(2) symmetry, our bounds 
severely constrain models of conformal technicolor. In = 1 superconformal theories, 
we place strong bounds on dim($^<l>), where <^ is a chiral operator. These bounds 
asymptote to the line dim(<I>1^<I>) < 2dim(<I>) near dim($) ~ 1, forbidding positive 
anomalous dimensions in this region. We also place novel upper and lower bounds on 
OPE coefficients of protected operators in the <I> x $ OPE. Finally, we find examples 
of lower bounds on central charges and flavor current two-point functions that scale 
with the size of global symmetry representations. In the case of M = 1 theories with 
an SU(iV) flavor symmetry, our bounds on current two-point functions lie within an 
0(1) factor of the values realized in super symmetric QCD in the conformal window. 
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1 Introduction 



Conformal phases in four dimensions are ubiquitous and may play a crucial role in beyond 
the Standard Model physics. Some simple examples include walking [1-6] or conformal [7-12] 
technicolor, explanations of the flavor hierarchies [13-16], and solutions to the supersymmet- 
ric flavor problem [17-30], the fi/Bfi problem in gauge mediation [31-36], or the rj problem in 
inflation [37]. Moreover, studying conformal field theories (CFTs) can also give us important 
insights into quantum gravity and string theory via the AdS/CFT correspondence [38-40], 
which in turn provides a simple framework for describing many new physics scenarios via 
effective field theories in AdS [41-43] (dual to 'effective CFTs' [44]). 

However, in recent years it has been realized that the restrictions imposed by conformal 
symmetry are not very well understood. While constraints on the form of simple correlation 
functions (e.g., [45, 46]) and unitarity restrictions on operator dimensions [47, 48] were 
worked out long ago, it was pointed out in [49] that crossing symmetry of four-point functions 
combined with the constraints of unitarity imply additional bounds on operator dimensions 
that must be satisfied in any consistent CFT. These bounds were soon strengthened [50] and 
extended to bounds on scalar operator product expansion (OPE) coefficients [51]. In [52] 
the bounds were also extended to A/" = 1 superconformal field theories (SCFTs); bounds 
on central charges in general CFTs and SCFTs were also explored in [52] and [53]. In 
addition, progress on incorporating global symmetries into the program (important for both 
phenomenological applications and to have a more direct comparison with known theories) 
was made in [54], and improved bounds (both for general CFTs with global symmetries and 
for SCFTs) were presented in [55]. 

The methods used in [49-55] to obtain bounds involve applying linear functionals to CFT 
crossing relations, which in practice means taking linear combinations of derivatives of the 
crossing relations evaluated at a particular point. By searching for linear functionals that are 
positive when acting on the contributions of all possible primary operators in the spectrum 
other than the unit operator, one can obtain bounds on OPE coefficients (and sometimes 
operator dimensions). However, to implement this positivity condition, the authors of [49-55] 
introduced a finely-spaced discretization of the set of possible operator dimensions, making 
the resulting linear programming problem numerically difficult and limiting how far the 
idea could be pushed. This numerical limitation was particularly apparent when considering 
systems of crossing relations that occur in theories with global symmetries, where the bounds 
obtained so far still seem to be quite far from their optimal values. 

In the present paper we will present an alternate approach that completely avoids this 
discretization of dimensions. We will use the fact that linear combinations of derivatives 
of conformal blocks can be arbitrarily-well approximated by ratios of polynomials in the 
operator dimensions, which allows us to convert the problem of obtaining bounds into a 
semidefinite programming problem that is numerically much more efficient. This then allows 
us to obtain much stronger bounds on CFTs and SCFTs, particularly in the presence of 
global symmetries. 

More concretely, for general CFTs we will consider four-point functions of scalar operators 
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0, as well as collections of operators (pi transforming as fundamentals under SO(A^) or SU(A^) 
global symmetries. For theories with A/" = 1 supersymmetry we will focus on the case of 
chiral superconformal primary operators $, as well as on collections of chiral operators $i 
transforming as SU(A^) fundamentals. We start by reviewing the relevant crossing relations 
and representation theory in section 2. There we will also introduce our new method to obtain 
bounds on operator dimensions and OPE coefficients based on semidefinite programming. 

In section 3 we use this method to derive general bounds on operator dimensions. In 
the case of general CFTs with SO(A^) global symmetries, we will place upper bounds on the 
dimension of the lowest-dimension SO(iV)-singlet operator appearing in the (pi^cpj OPE. This 
greatly improves upon the bounds in the presence of global symmetries previously presented 
in [54, 55]. We also place similar bounds on the lowest-dimension SO(A^) symmetric tensor 
(I>(i4>j)- In the case of SU(A^) global symmetries we can additionally place bounds on SU(A^)- 
singlet or SU(A^)-adjoint operators appearing in the 0j x 0^"'' OPE. Somewhat surprisingly, 
we find that SU(A^)-singlet bounds turn out to be identical to S0(2A^)-singlet bounds using 
the present method. 

The special case of an SO (4) or SU(2) global symmetry is relevant for the scenario of 
conformal technicolor [7], with or without custodial symmetry. In this scenario one would 
like the dimension of the Higgs operator H to be somewhat close to 1, while the dimension 
of H'^H should be close to or greater than 4. On the other hand, the bounds in this paper 
show that requiring dim(i7^i7) > 4 forces one to have at least dim(if) > 1.52, excluding 
fiavor-generic versions of this scenario and placing significant constraints on models where 
Yukawa-like suppressions are generated in four-fermion operators. 

In A/" = 1 superconformal theories we also place bounds on the lowest- dimension scalar 
superconformal primary appearing in the $ x $^ OPE, where $ is a chiral operator. This 
greatly strengthens the bounds presented in [52, 55]. In fact, we will see that the bound 
appears to asymptote to the line dim($"''$) < 2dim($) near dim($) ~ 1, essentially 
excluding the possibility of 'positive anomalous dimensions' (as recently discussed in [56]) 
in this region. This also implies that the solution to the /i/i?/i-problem proposed in [31, 32] 
cannot easily work near dim($) ~ 1. 

In section 4 we explore bounds on OPE coefficients. First we strengthen the upper 
bounds presented in [51] on the sizes of OPE coefficients of scalars O appearing in the 
(j) X (j) OPE in non-supersymmetric theories. Then, as a new application of these methods 
in superconformal theories, we place both upper and lower bounds the OPE coefficient of 
the chiral operator which always appears in the $ x $ OPE. In this case, lower bounds 
are possible because unitarity requires that there is a gap in the spectrum of dimensions, so 
no other nearby operators can mimic the effects of the $^ operator in the conformal block 
decomposition. We similarly place upper and lower bounds on the OPE coefficients of the 
other higher-spin protected operators that can appear in the $ x $ OPE. These bounds have 
interesting implications for Banks-Zaks theories or CFTs with weakly-coupled AdSs duals, 
where they can be checked in perturbation theory. 

Next, in section 5 we place lower bounds on the central charge c, which appears as the 
coefficient in the two-point function of the stress tensor: (TT) oc c. These bounds strengthen 
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and expand upon those previously explored in [52, 53, 55]. In theories with operators of 
dimension d transforming as fundamentals under SO(A^) or SU(A^) global symmetries, we 
find that the bounds scale linearly with N near ~ 1, consistent with our intuition from 
free CFTs. We explore these bounds on c in both general CFTs and Af = 1 SCFTs. In the 
latter case one can calculate c using 't Hooft anomaly matching in many known SCFTs, and 
our bounds are satisfied in all such examples that we have checked. 

In section 6 we place similar bounds on the coefficient k appearing in the two-point 
function of a global symmetry current: {J^J^) oc kTt{T^T^). Here we extend the previous 
results of [52] to include the full information about global symmetries. In the case of scalar 
operators transforming as fundamentals of SO(A^), we place lower bounds on kso{n)- In 
the case of SU(A^) global symmetries, one can either bound the OPE coefficient appearing 
in front of the SU(A^) (adjoint) current or the coefficient in front of an SU(A^)-singlet 
current corresponding to a different global symmetry. In the latter case, the bounds again 
scale linearly with near ~ 1 in accordance with our intuition from free CFTs. We 
also compute similar bounds in A/" = 1 SCFTs where k can be computed using 't Hooft 
anomaly matching, and present a comparison of our results with supersymmetric QCD in 
the conformal window [57]. We conclude in section 7. 



2 Bounds from Crossing Relations 
2.1 CFT Review 

Let us begin by reviewing some basic aspects of conformal field theories that will be important 
for our discussion. The conformal algebra contains, in addition to Poincare generators, a 
dilatation generator D and special conformal generators K^. Operators in a CFT can be 
classified into primaries satisfying Kf,O^{0) = 0, and their descendants P'' ■ ■ ■ P''O^{0)} 
Here, / denotes possible Lorentz indices. We will be primarily concerned with spin-£ 
operators which transform as traceless symmetric tensors of the Lorentz group, = O^^'"^'. 

Correlation functions of a conformal field theory on are completely determined by 
some simple discrete data: the spectrum of operator dimensions and spins, and the coeffi- 
cients appearing in the operator product expansion (OPE). Knowledge of the spectrum is 
sufficient to determine all two-point functions. For primary operators (9/ and with equal 
dimensions and spins {A,£}, we have 

{Oi{x^)0^{x,)) (X (2.1) 

where w^'^{x) is a tensor whose form is fixed by conformal symmetry (e.g., for spin-1 operators 
w'^^{x) = Tj'^^ — ^^^). When the dimensions and spins are not equal, the two-point function 

^Wc leave the adjoint action of charges on operators imphcit, i.e. K^O = [K^,0]. 
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must vanish. In addition, unitarity constrains A to satisfy [47, 48] 

A > 1 (^ = 0), 

A > £ + 2 (£>1). (2.2) 

These bounds can sometimes be strengthened if the conformal algebra is enhanced, as in 
superconformal theories. We will see some examples of this shortly. 

Let us choose an orthonormal basis of primaries Oi, so that the constant of proportionality 
in Eq. (2.1) is 6ij. Having done so, the remaining n-point functions of the theory are 
determined by coefficients in the operator product expansion. For real scalars (pi and 02, 
this takes the form [58] 

0i(x)02(O) = Yl h.^.oCiix,P)O\0), (2.3) 

where X,f,-i^,f,20 are constants that must be real in a unitary theory. The notation O ^ (pi x 
(p2 indicates that (9 is a primary operator in the OPE of (pi and (p2- We have grouped 
together each primary O and its descendants PO, P^O, . . . into a single term using the 
operator Ci{x, P) (which depends on the dimensions and spins of (pi, (p2, and O, though we 
are suppressing that dependence for brevity). One can show that the form of Ci{x,P) is 
completely fixed by conformal symmetry. For instance, applying special conformal generators 
Kfj^ to both sides of Eq. (2.3) gives a recursion relation for the terms in Cj which can be 
solved order- by-order. When (pi = (p2 = <P, Bose symmetry dictates that only even-spin 
operators may enter the OPE (2.3). 

In general field theories, the OPE is an asymptotic expansion, valid only at short dis- 
tances. However in a GET, because of the absence of scales, the OPE is an exact equality that 
can be used to simplify products of operators with arbitrary separation inside correlation 
functions.^ A key example for us is a four-point function of a scalar operator (p of dimension 
d, which can be evaluated as follows: 



(0(X1)</)(X2)0(X3)0(X4)) 



where we have inserted Eq. (2.3) twice and used Eq. (2.1) together with orthonormality of 
the O's. Here, A and i are the dimension and spin of O, u = \^ and v = l'^ f are con- 

^13^24 ^13^24 

formal cross-ratios, and the functions gA,eiu,v) are called conformal blocks. Since conformal 
symmetry completely fixes Ci and w^"^ , it also determines g^/. An exact expression in four 

^This is true provided there are no other operators 'nearby' in a sense that can be made precise. 



J2 \o\o'Ci{xi2,d2)Cj{xu,d,){0'{x2)0"{x,)) 



O'asxs 



Yl ^0-2lZ2d9AA^,v) 
■^12-^34 



Oecf>x 



xiixiiCi{xi2,d2)Cj{xs,,d,: 



w (2:24) 



A ' 

24 



(2.4) 
(2.5) 
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dimensions, computed by Dolan and Osborn [59, 60], is given by 



zz 

QaA'^^'") = -^^{kA+e{z)k^^i_2(z) - {z -H-z)), (2.6) 

z z 

kfsix) = x'^/2^Fi(/3/2,/3/2,/3,x), (2.7) 

where u = zz and v = [1 — z){l — ^) The unit operator is an important special case, with 
go,oiu,v) = 1. 



2.2 Crossing Relations for Singlets, SO(iV), and SU(iV) 

While a set of dimensions, spins, and OPE coefficients is enough to compute any correlation 
function, this data must satisfy additional consistency relations in a sensible CFT. To simplify 
(0000) using the OPE, we had to choose some way of pairing up the operators, and this 
choice necessarily broke manifest permutation symmetry among the 0(xj)'s. Nevertheless it 
should be the case that the end result remains permutation-symmetric, a requirement known 
as crossing symmetry. As an example, switching xi -f-)- in the conformal block expansion 
Eq. (2.5) leads to the crossing relation 

d 



^ X^gA/i'^^'^) = (^) ^o^A/(t^,M)- (2.8) 



0£<t>x<i> oe<f>x<i> 

Meanwhile, switching xi O X2 reproduces the statement that only even-spin primaries 
appear in x 0. Other permutations give no new information in this case.^ 

Recall that the Xo are real by unitarity, which means that the coefficients are 
nonnegative. This is a source of tension in Eq. (2.8), which can be expressed most clearly 
by rewriting our crossing relation as a 'sum rule' with positive coefficients, 

Fo,oiu,v)+ J2 ^o^A,£(«,^) = 0, (2.9) 

where 

FaAu,v) = , (2.10) 

Fo,oiu,v) = -1, (2.11) 

and we are suppressing the ^-dependence of F^/ for brevity. Note that we have isolated the 
term corresponding to the unit operator, whose OPE coefficient is fixed by the fact that has 
a canonically normalized two-point function. The unit operator contributes to Eq. (2.9) in 
its own particular way. Requiring that this contribution be cancelled by F^/s with positive 
coefficients leads to nontrivial constraints on the allowed A, i appearing in x 0. For some 



■^Our convention for conformal blocks here differs by a factor of (—2)^ from the one used in [59, 60]. 
''^Note that crossing symmetry of all four-point functions is equivalent to associativity of the OPE, which 
is enough to guarantee that higher n-point functions are crossing-symmetric as well. 
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explicit examples and many details about the structure of the sum rule, see [49, 50]. In 
sections 2.4 and 2.5, we will explain our improved method for extracting bounds on CFT 
data from Eq. (2.9). For now, let us present some generalizations of the sum rule for other 
kinds of operators. 



2.2.1 SO(iV) Crossing Relations 

An analysis of crossing relations in theories with SO(A^) and SU(A^) global symmetries was 
performed in [54], and improved bounds for SO(A^) were presented in [55]. We will make 
extensive use of these results, so let us review them here. 

Consider a real scalar primary 0j transforming in the fundamental representation of an 
SO(A^) global symmetry group. A complex scalar is a special case with symmetry group 
S0(2) = f/(l). Operators in 0j x 0^ can be organized into singlets S, symmetric tensors T, 
and antisymmetric tensors A of SO(A^). Schematically, 

0,x0, ~ ^5^^.(9 + + (2.12) 

S+ T+ A- 

The notation 5*^, T^, indicates that the sum is restricted to even-spin (+) or odd-spin (— ) 
primaries in (pi x (pj with the given representation, as dictated by Bose symmetry. Keeping 
track of the SO(A^) indices, each representation contributes differently to the conformal block 
decomposition of a four-point function. 



xl2xli{(pi{xi)(pj{x2)(t)k{,X'i)<pi{Xi)) = '^\oiSijSkl)gAA'^^'") 



s+ 

+ 5Z -^o {^ii^jk - SiJji) S'a.Km, v). (2.13) 

A- 

If we recompute this four-point function using a different operator pairing, each primary 
contributes again, but with the conformal cross-ratios u and v switched, and the tensor 
structures Sij6ki, Si^Sji, SaSjk permuted. Picking out the coefficients of each tensor structure 
then leads to three sum rules, which we can write in vectorial form 




0. (2.14) 



Here H/^^(^{u,v) is a symmetrized version of F/^^i{u,v), 



HaA^.v) = ^^j-^^ , (2.15) 

Ho,o{u,v) = 1. (2.16) 
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For brevity, we have not isolated the unit operator in Eq. (2.14); it is included with the 
even-spin singlets S~^. 

The case of SO (4) is special, since one can additionally decompose antisymmetric tensors 
into self-dual and anti-self-dual parts A±. Let us quickly summarize the consequences, 
though they will turn out to be irrelevant for this work. A new tensor structure can now 
appear in , namely eyki- In Eq. (2.13) we must replace 

^Xa{Sii6jk - SikSji) gA,e{u,v) ^XoiSuSjk - SikSji±eijki) qaA^^^)- (2-17) 

Since eijki maps to itself under permutations, the sum rule Eq. (2.14) is unaffected. We must 
simply supplement it with 

J2xlFA,i-J2^oFA,i = 0. (2.18) 



Before we proceed, it is worth mentioning that all three of the sum rules given in Eq. (2. 14) 
can be derived from a single 'master' crossing relation 



S+ T+ \ T+ A- 



Adding Eq. (2.19) to itself with u v gives the second row of Eq. (2.14) and subtracting it 
from itself with u ^ v gives the third row. To obtain the first row, we must make repeated 
use of the identity g{u,v) = {—lYg{u/v, 1/f): 



A- 



T+ A- T+ 



.d 



\ S+ T+ 



vJ 

d 



^gAA^/^^ + Y9aA^/^^ 

T+ A- 



T+ A- 

This implies in particular that the first sum rule is not independent from the other two. 
However, in practice we find it useful to retain all three sum rules, since we will keep only 
a finite number of terms in their Taylor expansions around a single point in {u, w)-space. 
Since the derivation Eq. (2.20) requires transformation between different {u,v) points, the 
exact equivalence between the third sum rule and the other two is only visible with an 
infinite number of terms in the Taylor expansion. However, it will be important to clarify 
the meaning of this 'master' sum rule (and its generalization to other symmetries) in future 
studies. 
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2.2.2 SU(iV) Crossing Relations 



Let us now consider a complex scalar 0^ transforming in the fundamental representation 
of an SU(A^) global symmetry. For this paper, we will only analyze four-point functions 
(0j(/)^Vfc0'^) that would be invariant under an additional U{1) acting on (p. Note that this is 
not tantamount to assuming such a U{1) exists — rather, we are restricting our attention to 
a subset of CFT correlators. The various channels for decomposing our four-point function 
now involve two different kinds of OPEs. Firstly, 



X 0- 



(2.21) 



which can contain SU(A^) singlets and adjoints of any spin. We also have 



(2.22) 



T+ 



containing symmetric and antisymmetric tensors with even and odd spins, respectively, and 
its complex conjugate 0*^ x 0^''' containing the conjugate operators in dual representations. 
Extracting the coefficients of different tensor structures in all possible ways of evaluating 



leads to the six-fold sum rule 



A- 



(2.23) 



T+ 



where 



V 



5± 



/ F \ 




H 




i-YF 
i-YH 









V J 





-(1 + ^)^ 

-i-y^F 

i-YF 
V i-YH I 



( \ 



F 

-H 
F 

\-H/ 



/ \ 


F 
-H 
-F 
\ H ) 



.(2.24) 



Once again, the unit operator is included among even-spin singlets . 



2.3 Crossing Relations in Superconformal Theories 

The 4D M = \ superconformal algebra extends the conformal algebra to include supersym- 
metry generators Qa, Qq, superconformal generators S'q, 5^, and a f/(l) -R-charge generator. 
SOFT operators admit a more refined classification into superconformal primaries satisfying 
SO{0) = SO{0) = 0, with their superconformal descendants obtained by acting with any 
combination of Q, Q, and P. It's easy to see using {S, S} ~ K that a superconformal primary 
is also a conformal primary. But the converse is not necessarily true. A multiplet built from 
a single superconformal primary generally contains several (though finitely many) conformal 
primaries whose dimensions, spins, and OPE coefficients are related by supersymmetry. 
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A principle example for this work is a cliiral superconformal primary scalar $ of dimen- 
sion which satisfies S'<l>(0) = S'$(0) = (5$(0). Unitarity implies that its dimension is 
proportional to its -R-charge, d = Below we will review the structure of the OPEs 

needed to decompose four-point functions of $ and $^ into conformal blocks. We also refer 
the reader to [52, 55, 61] for additional discussions of these OPEs. 

First, superconformal primaries O appearing in $ x are restricted to have vanishing 
i?-charge and a dimension satisfying the unitarity bound A > 2 + i, where £ > is 
the spin of O. Each superconformal primary generically comes with three superconformal 
descendants of definite spin which are also primaries under the conformal subalgebra. (When 
the unitarity bound is saturated, A = 2+i, two of these descendants vanish and the multiplet 
is shortened.) Schematically, the OPE takes the form 

<|.x<l>t ~ J2 [o + {QQO)e_i + {QQO),+i + Q^Q^oj, (2.25) 

where (9 G $ x denotes that the sum is over superconfoTmal primaries in $ x and the 
subscript on QQO indicates the spin. We are being somewhat sketchy in our notation; 
the exact form of these conformal primaries depends on A and i and is given in [52]. 
Superconformal symmetry imposes the following relations between their OPE coefficients,^ 

^{QQOWi = 4(A+^ + l)^^' ^^-^^^ 

^((3QO),_i = i(A-£-l)^^' ^^'^^"^ 
- (A + ^)(A-£-2) 
\'Q'o - 16(A + £+l)(A-^-l)^^- ^^-^^^ 

Note that A^„7=r^^, and „—2 vanish when A = i + 2, consistent with shortening of the 
superconformal multiplet. 

Meanwhile, the $ x $ OPE can only contain operators which are killed by Q. First 
and foremost, we have the chiral primary whose dimension is exactly 2d, by virtue of 
the relation between dimension and -R-charge for chiral operators. All other operators are 
Q-descendants. Schematically, 



$x<l> ~ $2+ QOe + ^QO. (2.29) 

i=2,4,... O 

The operators Oe transform in (|, ■^) representations of the Lorentz group S0(4) = SU(2) x 
SU(2), and satisfy the BPS shortening condition Q^(9^"3---«<',«i - af _ rjj^^ product QOi is 
then a spin-£ operator, which is required by the superconformal algebra to have dimension 

2 

2d + i. Finally, the remaining operators Q O are not protected by a BPS condition, and 
can have any dimension satisfying A > \2d — ?)\ + ?> + ^. Note that when d < 3/2, a 



^The difference in normalization from the formulae in [52] is due to our different convention for conformal 
blocks Eq. (2.6). 
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gap in dimensions exists between the protected operators ^^,QOe and the non-protected 
operators Q O. In contrast to the situation for $ x each conformal primary in Eq. (2.29) 
appears with an independent coefficient — there are no additional relations imposed by 
supersymmetry among operators in $ x 

Because of the U{\)r symmetry, crossing symmetry of the four-point function ($$"''$$"'') 
is a special case of crossing symmetry for SO (2). Note that the antisymmetric tensor 
representation of SO (2) is the trivial representation, so that we may equivalently write S~ 
(odd-spin singlets) for A~ (odd-spin antisymmetric tensors). We are also free to multiply 
the sum rule by any invertible matrix without changing its content.^ Consequently, we can 
rewrite Eq. (2.14) for SO (2) as 



5± V (-)'^A/ I \ -HA,i J 

In our superconformal four-point function ($$'''$<I>^^), the 5*^ terms will come from the OPE 
(2.25), while the T"*" terms come from (2.29). Making use of the relations between OPE 
coefficients Eq. (2.26), the specialization of Eq. (2.30) to the superconformal case is 




0, (2.31) 



where 



_ ^ (A + i) (A -£-2) 

J'a, = Fa, + + ^ + ^f^^W + 4(A - £ - 1)^^+^'^"^ 

'''''' ^ -Fa+2,1- (2.32) 



16(A + £+l)(A-£-l)' 



In addition, is J-" with odd spins flipped Fa/ — )■ {—YFa/ throughout, and Ti is with 
Fa/ -> Ha/- The set Tgpg consists of the BPS operators appearing in $ x $, namely <I>^ 
and QOe for i E {2,4,...}. Tj^j^_gpg consists of the remaining operators in $ x $. In going 
from Eq. (2.30) to Eq. (2.31), we have removed the factors of 2 in front of the symmetric 
tensor contributions because their conventional normalization differs between SO (2) and 
[^T+]u{i) = 2[A^+]sO(2)- 



2.3.1 Superconformal SU(A^) Crossing Relations 

It is straightforward to generalize this analysis to the case of a scalar superconformal primary 
$j transforming as a fundamental under an SU(A^) global symmetry. The index structure of 
the OPE is the same as is given in Eqs. (2.21) and (2.22), but with the additional constraints 

^Specifically, we will replace the middle row with itself plus twice the top row, and the top row with the 
middle row. 
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imposed by supersymmetry discussed above. Note that now both BPS and non-BPS odd- 
spin operators can appear in the $j x $j OPE as SU(A^) antisymmetric tensors. Including 
these constraints, the six-fold sum rule of Eq. (2.23) becomes 



S± Ad± 



A- 

2d+e,£ 



+ 



)i2 T/ 



T+ 



+ 



(2.33) 
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(2.34) 



2.4 Bounds from Crossing Relations 

Crossing symmetry of four-point functions encodes an infinite number of relations between 
OPE coefficients — one for each value of the conformal cross-ratios u and v. In [49] a general 
method was outhned for extracting bounds on CFT data using these relations, together with 
the constraints of unitarity. We will now review this method for the simplest case of a real 
scalar of dimension d. Subsequently, we will discuss how the original method can be 
improved using semidefinite programming. 

Suppose we would like to bound the OPE coefficient of a particular operator Oq of 
dimension Aq and spin appearing in </> x 0. The first step is to isolate A^^^ on one side of 
the sum rule Eq. (2.9), 

Ao„Fao,£o(^>^) = -i^o,o(^^,f) - ^ A0FA/(n,w). (2.35) 

We can obtain different expressions for A^^^ in terms of the other OPE coefficients by 
evaluating Eq. (2.35) at different values of u and v. We could also take some number of 
u- and f-derivatives first, and then evaluate. And in general, we can apply any linear 
functional a to both sides, 

A^,a(FA„,,o) = -a(Fo,o) - A^a(FA,,). (2.36) 

A key insight of [49] is that the functions Fa/ share certain positivity properties, so that it's 
sometimes possible to find a linear functional a such that 

a(-^Ao/o) = and (2.37) 

o(-Pa,^) > 0, for all other (non-unit) operators in the spectrum. (2.38) 
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Eq. (2.37) is simply a normalization condition, but to satisfy Eq. (2.38) one must choose 
a carefully. If a satisfies these constraints, then since the are positive by unitarity, 
Eq. (2.36) becomes an upper bound on A^^^, 

A^^ = -a(Fo,o) - P°^- ^ P°^- ^ -"(^o,o)- (2.39) 

The space of viable a's depends on precisely what assumptions one makes about the spectrum 
of the CFT. If one makes an assumption about the spectrum of operator dimensions that 
makes it easier to satisfy Eq. (2.38) (e.g., all scalars have a dimension greater than some 
Amin) and then finds a linear functional a such that the bound of Eq. (2.39) violates the 
unitarity constraint A^^^ > 0, one can rule out that assumption about the spectrum. 

Now, to make the bound (2.39) as strong as possible, we should minimize — a(Fo^o) over 
the set S of all a satisfying the constraints (2.37, 2.38). These constraints carve out a 
convex subset of the space of linear functionals, so the task of determining the best a is 
an infinite-dimensional convex optimization problem. It would be extremely interesting to 
develop analytical techniques for finding solutions. However, the most successful approaches 
to date, including the one we present here, involve simplifying the problem to make it 
tractable on a computer, and then determining solutions numerically. 

Putting our optimization problem on a computer requires surmounting two difficulties: 

1. The search space S of a's satisfying Eqs. (2.37, 2.38) is infinite dimensional. 

2. The number of constraints a{FA,e) > is infinite — there's one for each A,i. 

The first difficulty is easy enough to address: we can restrict to a finite-dimensional 
subspace W of linear functionals. Then, minimizing — Q!(Fo,o) over all a G W fl 5 will give 
a possibly sub-optimal, but still valid bound A^^^ < — a(Fo^o)- The choice of W is somewhat 
arbitrary, and it would be interesting to explore a wider variety of functionals than we do 
here. Following [49-55], we will simply take linear combinations of derivatives around the 
symmetric point z = 1 = 1/2. That is, we define Wk to be the space of functionals 

a:Fiz,z) ^ a^„ar9^F(l/2, 1/2), (2.40) 

m+n<2k 

with real coefficients amn- This choice is computationally convenient, and will prove useful 
in our solution to the second difficulty in a moment. One hopes that as we increase k to 
include more and more derivatives, our search will cover more and more of S, and our bound 
will converge to the optimal one. 

The second difficulty is more problematic. Since angular momentum i is discrete, it's 
reasonable to include constraints with i = 0,1, . . . , L, for some large L. But the dimension 
A can vary continuously, and the constraints q;(Fa,^) > carve out a complicated shape S 
inside W as A varies. The computer has to know about this shape, which means we must 
encode it with some finite amount of data. The approach used in [49-55] is to approximate 
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«(^Ao/o) = 1 



"(^Aoa) = 1 



a{FA,,i,) > 



Figure 1: The 'search space' 5 (shown in blue) is the intersection of the hyperplane a(-FAo,^o) ~ ^ 
with the convex cone of hnear functionals a satisfying a(-FA,£) ^ for all {A,£} in the spectrum. 
Previous methods discretized A to some finite set {Aj}, thus approximating S as an intersection 
of a finite number of hyperplanes and half-spaces (left). Our approach is to approximate S as 
the intersection of a smaller number of curved spaces — specifically cones of semidefinite matrices 
(right). Such intersections are sometimes called spectrahedra. 

the shape by a convex polytope — namely discretize A to lie in some finite set {Aj}, so 
that the constraints a{F/^.^£) > become a finite number of linear inequalities for a. Then 
the problem of minimizing — a(Fo,o) becomes a linear programming problem, which can be 
solved by jumping from vertex to vertex on the boundary of the polytope, following the 
direction of steepest descent. As one makes the set {Aj} larger, the approximation of S as 
a polytope gets more and more refined, and the solution should converge to the correct one. 

This method can be quite powerful if one chooses the {Aj} carefully. However, some 
basic tensions limit how far it can be pushed. For example, consider increasing k to obtain 
a stronger bound. At higher k, the space Wfc can include wilder linear functionals, and 
one must include more Aj to ensure that a constraint <y{F/^j) > isn't violated. However, 
the running time of the usual search algorithm is cubic in the number of constraints, which 
means that computations become quickly unwieldy. 

Our approach in the present paper is to approximate S with a different kind of shape 
that is more efficient to encode than a polytope, one that naturally respects the properties 
of conformal blocks (specifically the differential equation that they satisfy), and also admits 
fast searches. In the process, we will do away with the discretization A G {Aj} entirely. 

2.5 Semidefinite Programming 

Semidefinite programs (SDPs) [62] are linear optimization problems that can contain positive- 
semidefiniteness constraints for matrices, along with the usual linear inequalities included 
in linear programs. As we'll see momentarily, positive-semidefiniteness lets us express the 
condition that a collection of polynomials be nonnegative for all values of their arguments. 
This is useful for us because there is a systematic approximation for the derivatives of Fa,z 
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in terms of polynomials. Specifically, there exist positive functions Xii'^) polynomials 
PP'"'{A) such that 

ara^FA,,(l/2, 1/2) ^ xKA)Pr(A), (2.41) 

where the approximation can be made arbitrarily good, at the cost of increasing the degree 
of Pp^. The details of this approximation, which follows from the differential equation for 
conformal blocks along with some basic facts about hypergeometric functions, are explained 
in appendix A. 

For now, let us assume that such an approximation exists, and understand how to phrase 
our problem as an SDP. We will write F/""'(A) = ^^S^Fa^^ (1/2, 1/2) for brevity. Once 
again, we would like to minimize —amnF^"'{0) subject to the constraints 

a„„F,7(Ao) = 1, (2.42) 
amn^r"(A) > for A > A^, for all < £ < L, (2.43) 

where A^ is a lower bound on A depending on the spin i. 

Using Eq. (2.41) along with the fact that Xii^) is positive, Eq. (2.43) becomes the 
statement that each polynomial amn-P/""'(Af (l+x)) is nonnegative on the interval x G [0, oo). 
Such statements are naturally written in terms of positive-semidefinite matrices, a fact which 
is well-known in the optimization literature and has been exploited to solve a wide variety 
of problems (see, e.g., [63]). The rewriting proceeds as follows. Firstly, a theorem due to 
Hilbert [64] states that a polynomial p{x) is nonnegative on [0, oo) if and only if 

p{x) = f{x) + xg{x), (2.44) 

where both f{x) and g{x) are sums of squares of polynomials. Now suppose / and g have 
degrees 2d and 2d' respectively, and let [x]d denote the vector with entries (l,x, . . . ,x'^). If 
f{x) is a sum of squares of polynomials with coefficients Cj = (qo, . . . , Q^), then we have 




(2.45) 



where A = - Cjcf is positive-semidefinite. Conversely, any positive-semidefinite matrix A 
admits a Cholesky decomposition A = Cjcf , so that [x]j74[x]d is a sum of squares. Thus, 
the condition that p{x) be nonnegative on [0, oo) can be written 

p{x) = [x]^A[xU + x{[x]lB[xU,), with A,BhO, (2.46) 

where the >z symbol means 'positive-semidefinite.' 

Returning to OPE bounds, we now have the following presentation of our convex opti- 
mization problem as an SDP: minimize —amnF^^{0), subject to the constraints 

(2.47) 

for < £ < L, (2.48) 
for < £ < L. (2.49) 



a^„Pr"(A,(l + x)) = [x]lA,[xU, + x{[x]lBe[xU'J 
Ai, BehO 
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There are numerous advantages to this formulation. Firstly, we avoid discretizing the set 
of operator dimensions A, and thus evade the trade-off between refining {Aj} and improving 
the running time. Further, small and large A are accounted for equally well, so there is no 
need for separate checks on the asymptotic behavior of a (-Fa/) at large dimensions. Most 
importantly, there exist efficient algorithms for solving semidefinite programs using interior 
point methods, with some excellent implementations (see appendix B). Their complexity 
scales much less sharply with the dimension of the search space than the linear programming 
algorithms used in [49-55]. Consequently, we have been able to push the previous state-of- 
the-art searches from 55 dimensions to almost 400 dimensions in some cases. 



2.6 Generalizations for Global Symmetries 

While we specialized the above discussion to the case of the singlet sum rule Eq. (2.9), it 
is straightforward to modify it for situations with global symmetries. E.g., if we wish to 
place a bound on the OPE coefficient of an operator appearing in the SO(A^) sum rule 
of Eq. (2.14), we should look for a vectorial linear functional a satisfying 

= 1, (2.50) 

> 0, for all other (non-unit) operators in , (2.51) 
a \ [l — j^) Fa/ I > 0, for all operators in T^, and (2.52) 

> 0, for all operators in A~ . (2.53) 




Any such linear functional then leads to the upper bound 





AL < -« I Fo,o I . (2.54) 



The modification for alternatively placing bounds on the OPE coefficients of or A~ 
operators should be clear. As in the singlet case, we can also rule out an assumption about 
the spectrum of operator dimensions by making the assumption and then finding a linear 
functional that leads to a violation of the unitarity constraint A^^^ > 0. 

Similarly, we can bound the OPE coefficient of an operator appearing in the SU(A^) 
sum rule of Eq. (2.23) by finding an a satisfying 



a 

a 



(Co) = 1> (2-55) 
i^A/) — 0' for all other (non-unit) operators in the spectrum, (2.56) 
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where / = {S^ , Ad^ ,T+ , A"}. This leads to the upper bound 



(2.57) 



The appropriate generahzation of this logic for placing bounds on operators in other SU(A^) 
representations, and also for obtaining bounds using the superconformal sum rules given in 
Eqs. (2.31) and (2.33), should be clear. 

In all of these situations, the task of numerically finding the optimal a can be recast in 
terms of a semidefinite program. Similar to what we described in the previous section, to do 
this we use the fact that derivatives of any of the functions {-Fa/, -f^A/, J^A,e, 'Ha/, J^a/, 'Ha/} 
at (1/2,1/2) can be arbitrarily-well approximated by positive functions times polynomials 
in A. The details of these approximations can be found in appendix A. 



2.7 Coincidence Between SU(iV) and S0(2iV) Singlet Bounds 



In the course of running the above algorithm, we found that our bounds on singlet operators 
appearing in an OPE between SU(A^) fundamentals were numerically identical to bounds 
on singlets appearing in an OPE between S0(2A^) fundamentals. This exact coincidence 
is surprising given the rather different structure of the crossing symmetry constraints. It 
hadn't been previously observed because SU(A^) computations were too difficult to perform 
with previous techniques. In this section, we'll discuss the relations between those bounds 
in more detail. 

Let us consider more generally a GET with global symmetry group Q. Suppose we 
want to obtain a dimension bound on a singlet scalar operator entering a given OPE. The Q 
crossing symmetry constraints produce a bound Ag{d). Now consider a subgroup H d Q and 
repeat the procedure. This time, the l-L crossing symmetry constraints will produce a bound 
Ay^{d). At this point we must distinguish two cases: 1) all H-singlets are also ^-singlets, 2) 
some nontrivial representation of Q, once decomposed with respect to the subgroup, contains 
"H-singlets. In the first case we can immediately conclude 



The above inequality is clear: there are no CFT's with global symmetry "H where the first 
scalar singlet operator entering a given OPE has dimension larger than A^((i). Thus in 
particular there are no CFT's with a larger global symmetry. 

An example of such a group and subgroup is given precisely by SU(A^) C S0(2A^). In the 
decomposition with respect to the subgroup, the only singlets come from S0(2A^)-singlets: 
the symmetric tensor goes to a symmetric tensor and an adjoint while the antisymmetric 
tensor goes to an antisymmetric tensor and an adjoint. Thus, it is natural to expect the 
triple sum rule Eq. (2.14) to give a bound stronger than or equal to the sextuple sum rule 
Eq. (2.23). Indeed, one can verify this explicitly at the level of the optimization problem for 
a. 



Ag{d)<An{d) 



(^-bound stronger). 



(2.58) 
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To prove the equality of SU(A^) and S0(2A^) bounds one should also show that whenever 
a linear functional satisfying Eqs. (2.50-2.53) exists, it is possible to construct a second linear 
functional satisfying Eq. (2.55). Unfortunately, we have not been able to find an analytic 
proof of this result. However, we find numerically that it is always possible - it would be 
good in future studies to gain a deeper understanding of why this is the case. 

In the case 2) the two bounds are unrelated, since the H-bound could in principle be 
determined by representations coming from the decomposition of nontrivial representations 
of the larger symmetry group. This is the case for SO(A^) and SO(A^') or SU(A^) and SU(A^'), 
with > A^'. In these examples we numerically observe behavior opposite to (2.58). 

3 Bounds on Operator Dimensions 

3.1 General Theories 

As a first application of our semidefinite programming algorithm, let us reproduce the singlet 
dimension bound first derived in [49], and later improved in [50]. We let be a real scalar 
of dimension in a general CFT, and seek to place an upper bound on the dimension of 
0^, the lowest dimension scalar appearing in x 0. The procedure is precisely as described 
in section 2.4. In figure 2, we show the resulting bounds for k = 2, ... ,11, with /c = 10 (a 
5 5- dimensional search-space) being the previous state-of-the-art. We find perfect agreement 
with older linear programming-based calculations for each k = 2, ... ,10. 

The curves appear to converge at large k, which is perhaps indicative that they are 
approaching the best possible bound given our assumptions (referred to as foo{,d) in [50]).^ 
We will see this kind of convergence in many other plots in this paper. An approximate fit 
to the strongest {k = 11) bound is given by^ 

dim(02) < 2 + 3.006e + 0.160(1 -e-^°'), (3.1) 

where d = 1 + e, with e between and 1. Notice that the behavior for both small and large 
e is approximately linear. The bound crosses dim(0^) = 4 around d ~ 1.61. 

3.2 Singlet Operators in SO(iV) and SU(iV) Theories 

We can also place bounds on the lowest dimension singlet appearing in 0j x 0^, where 0, 
transforms as a vector of an SO(A^) global symmetry. The procedure is as described in 
section 2.6, where we must assume that A > Amin for all scalars in , and then scan over 
Amin to obtain a dimension bound. Recall from section 2.7 that our bounds on singlets of 

^However, since the full optimization problem involves an infinite-dimensional search space, it's always 
possible a new search direction could open up at higher k. Fully establishing convergence would require more 
detailed analysis than we do here. 

®While it gives a good description of the shape, we have chosen this functional form somewhat arbitrarily; 
it is possible that a different basis of functions should be used when describing the optimal bound. 
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Upper bound on dim(<^^) 




Figure 2: An upper bound on the dimension of (/>^, the lowest dimension scalar appearing in (j) x cj). 
Curves for k = 2, ... ,11 are shown, with the k = 11 bound being the strongest. 



SU(A^) turn out to be identical to those for singlets of S0(2A^). Hence, we will present all 
SU and SO singlet bounds together, with even values of standing for both SO(A^) and 
SU(iV/2). 

Previous attempts to compute bounds for theories with global symmetries have been 
somewhat hindered by the need to optimize over very high-dimensional spaces. Since the 
vectorial sum rule Eq. (2.14) has three components, a given k corresponds to 

k(k + l) , , 

X 3 (3.2) 

different linear functionals. The linear programming methods implemented so far are essen- 
tially limited to a search space dimension that is not much larger than ~ 50, or /c ~ 5 for 
SO(A^). Worse, SU(A^) vectorial sum rules have six components, making them even harder 
to explore. However, our semidefinite programming algorithm appears to have few problems 
with large search spaces, and we will present most of our bounds up to A; = 11, regardless of 
the type of global symmetry group. 

As an example, figure 3 shows a bound on the lowest dimension singlet in theories with 
an SU(2) or S0(4) global symmetry.^ This bound is particularly interesting for conformal 

^Note that to compute the S0(4) bound, we have only used the triple sum rule of Eq. (2.14). It 
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technicolor models, as we will discuss in detail in the following section. Notice again that 
the curves start to converge at large k. An approximate fit to the strongest {k = 11) bound 
is given by 

dim{\(f)\^) < 2 + 3.119e + 0.398(l (3.3) 
where d = 1 + e, with e between and 1. This bound crosses Aq = 4 around d ~ 1.52. 



Upper bound on dim{\(j)\'^) for S0(4) or SU(2) 

Ao 




Figure 3: An upper bound on the dimension of (p^cj), the lowest dimension singlet scalar appearing 
in (j)^ X (j), where (f) transforms in the fundamental representation of an S0(4) or an SU(2) global 
symmetry. Curves are shown for A; = 2, . . . , 11. The bounds for S0(4) and SU(2) are identical in 
each case. The strongest bound crosses Aq = 4 around d = 1.52. 

Figure 4 shows dimension bounds for SO(A^) with = 2, . . . , 14 and SU(A^) with = 
2, . . . , 7. The strongest bound corresponds to the global symmetry group S0(2) = U{1), 
and the bounds weaken as A^ increases. One might naively expect a larger symmetry group 
to produce a stronger bound. For instance, a theory with an SO(A^) symmetry certainly 
also has an SO(A^ — 1) symmetry, so why shouldn't all bounds from the former apply to the 
latter? However, as discussed in section 2.7, the problem we are solving actually changes 
with A^, and this turns out to be a more important effect than the enhanced symmetry. Note 
that the lowest dimension singlet under an SO(A^ — 1) subgroup of SO(A^) is not necessarily a 
singlet at ah under the full SO(A^). Thus, SO(A^) bounds for larger A^ apply to the operator 

is straightforward to verify that including the fourth sum rule of Eq. (2.18) leads to a redundant set of 
constraints, and is therefore unnecessary. 
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with lowest dimension among a more restricted class of operators, and consequently can be 
weaker. 



Upper bound on dim(|0|2) for SO(iV) or SU(iV/2), = 2, . . . , 14 

Ao 




Figure 4: An upper bound on the dimension of the lowest dimension singlet scalar appearing 
in (fy'" X (j) (or (j)^ x (/)), where (j) transforms in the fundamental representation of an SO(iV) global 
symmetry or an SU(iV/2) global symmetry (when iV > 4 is even). Curves are shown for N = 
2, . . . , 14, with N = 2 being the strongest bound. 



3.2.1 Implications for Conformal Technicolor 

Let us briefly discuss some phenomenological implications of the bounds presented in figures 3 
and 4. A more detailed discussion of these implications will also appear in [65], and our 
analysis draws heavily on the previous discussions of [7-10, 49, 54, 55], as well as the recent 
talk of [66]. 

Arguably the most interesting operator dimension in the Standard Model is dim(if^if), 
the dimension of the Higgs mass operator, where H transforms as a bifundamental under 
SU(2)i X U(1)y C SU(2)l X 811(2)/?. In a weakly-coupled theory with a scalar Higgs, 
this dimension is approximately 2, which leads to the hierarchy problem and its associated 
puzzles. 

The idea of increasing dim(if^^i/) to ameliorate the hierarchy problem is an old one. 
In traditional Technicolor models, the role of the Higgs is played by a fermion condensate 
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ipip with dimension 3, so that the 'mass' term {ipipY is irrelevant. A basic tension in this 
setup is that the 'Yukawa' terms {ipip)qu which generate fermion masses after EWSB are 
also irrelevant. To correctly account for the top-mass, we must imagine that such terms 
are suppressed by a low scale in the Lagrangian £Yuk. 3 a^(^^)^^- -^^^ ^^^^ same low 

low 

scale would then generically appear in other four-fermion operators, leading to dangerous 
flavor- changing neutral currents. 

Conformal Technicolor (CTC) [7] seeks to avoid this tension by assuming that H par- 
ticipates in strong conformal dynamics above the electroweak scale, which generates a large 
dimension for H'^H, while the dimension of H remains near 1. While this idea is intriguing, 
we will show that it needs additional assumptions to work in practice. In particular, our 
bounds definitively rule out the simplest 'flavor-generic' CTC models. 

To begin, let us determine the range of d = dim{H) and A = dim{H^H) that is 
phenomenologically viable in CTC. Firstly, we must require that yt remain perturbative 
throughout the conformal regime, which places an upper bound on the possible running 
distance. Indeed, suppose conformal dynamics occurs between Aew ~ 47rw and some higher 
scale Auv- Within this range of energies, Yukawa couplings run according to 



mif^) = hr- ?/^(^Ew) (3.4) 

(ignoring corrections from small perturbations away from exact conformal symmetry, like 
SM gauge couplings and other Yukawa couplings). Requiring yt < 47r for all /i G [Aew, Auv] 
then gives 



1 

Auv ^ I ^Ew\''-^ 



Aew ^ \ rrit J 



(3.5) 



Secondly, we must ensure that small perturbations of the theory by the Higgs mass 
operator H'^H don't destabilize the conformal dynamics. This is certainly the case if WH 
is irrelevant, A > 4. If on the other hand A < 4, then we must also impose the lower bound. 



Auv < ^ T 



Aew Vc(A 



(3.6) 



where c(Auv) is the coefficient of in the perturbation 5C = c{A\jy)WH at Auv- The 
strength of the bound Eq. (3.6) varies, depending on the amount of tuning we're willing to 
tolerate in this coefficient. 

Finally, while Eqs. (3.5) and (3.6) prefer a small running distance, Auv must also be 
sufficiently large to suppress problematic flavor-changing operators, such as {ds'^) (sd'^) which 
contributes to K-K mixing. In a 'flavor-generic' model, we should demand 

1 1 

~ TT (generically), (3.7) 



A2 ~ A2 
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where Ap ~ 3.2 x lO^TeV for CP-violating contributions to K-K mixing [67]. More 
optimistically, we might imagine that {ds'^){sd'^) is generated with Yukawa suppression, so 
that the constraint above gets modified to 



yd(Auv)2/^(Auv) 




(optimistically) , 



(3.8) 



with ?/j(Auv) given by Eq. (3.4). 

Together, these requirements restrict viable models to a particular region of the d-A plane, 
which can then be compared with our bounds. In models where the conformal dynamics is 
custodially-symmetric, H transforms as a fundamental of S0(4) = SU(2)i x SU(2)ii; (which 
is weakly gauged by SM gauge fields). However, the assumption of custodial symmetry is 
not actually necessary for us because our bound for SU(2)i alone is identical to our bound 



The viable regions for fiavor-generic and fiavor-optimistic CTC models are shown in 
figure 5, superimposed with our strongest SU(2) dimension bound. The right-hand edge of 
the viable regions comes from the combination of Eq. (3.7) with Eq. (3.5), while the bottom 
edges come from the combination of Eq. (3.7) with Eq. (3.6) for different values of c(Auv). 
We see that for reasonable assumptions about the coefficient c(Auv), fiavor-generic models 
are ruled out. This conclusion remains true even if the conformal dynamics respects CP 
symmetry, in which case the effective flavor scale can be closer to Air ~ 10^ TeV. 

By contrast, flavor-optimistic models with reasonable tunings c(Auv) ^ 0.1 and some- 
what large dimensions d ~ 1.3-1.5 are not necessarily ruled out. Our bound does place an 
upper limit on the scale of new physics Auv, but with sufficient Yukawa suppression these 
upper limits can be phenomenologically acceptable. For instance, with c = 0.01, Auv niust 
lie below 6.8 x 10'^ TeV, while c = 0.1 gives Auv ^ 1-6 x 10'^ TeV. At some point however, 
the predictions for these models become essentially those of minimal fiavor violation with 
a low fiavor scale, and strong conformal dynamics seems more and more like a gratuitous 
assumption. 

3.3 Symmetric Tensors in SO(A^) Theories 

It is straightforward to modify our procedure to obtain bounds on symmetric tensors 
appearing in 0j x (pj. To bound a symmetric tensor with dimension Aq and spin Iq, we look 
for a linear functional satisfying the normalization condition 



as well as <y{V) > for all other vectors V in the SO(A^) sum rule. 

Figure 6 shows the resulting dimension bound on (the lowest dimension scalar 

symmetric tensor appearing in 0j x (pj) in the case of S0(4) symmetry. Note that this bound 



for SO (4). 




1, 



(3.9) 
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Viable regions for Conformal Technicolor models 




Figure 5: Viable regions for conformal technicolor models in the flavor-generic (red) and flavor- 
optimistic (cross-hatched green) cases are shown superimposed with our bound (blue, excluding the 
gray-shaded region). Regions for c(Auv) = 1) 0.1, and 0.01 are shown in successively lighter shades 
of each color, with the largest region corresponding to c(Auv) = 0.01 in each case. Flavor-generic 
models are ruled out. 

does not apply in a simple way to operators in theories with SU(2) symmetries, because 
there is no coincidence between SU(A^) and S0(2A^) bounds for non-singlets. 

3.4 Superconformal Theories 

Now let us turn to bounding operator dimensions in superconformal theories, using the sum 
rule Eq. (2.31). A bound on dim($"''$) in terms of dim($) was first obtained in [52] using 
only the middle row of Eq. (2.31). In [55], it was shown that the bound could be improved by 
incorporating the other rows, and linear programming calculations were given up to /c = 4. 
In figure 7, we present a new version of these bounds for k up to 11, corresponding to a 
19 8- dimensional search space. 

Several interesting new features emerge at large k. Most strikingly, the bound appears to 
be tangent to the factorization line Aq = 2d near d = 1. Figure 8 shows a higher- resolution 
plot for small values of d, which displays this behavior more clearly. An approximate fit to 
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Upper bound on dim(0(j0j)) for SO (4) 




Figure 6: An upper bound on the lowest dimension symmetric tensor scalar appearing in (j) x (j), 
where (p transforms in the fundamental of SO (4). Here we show k = 2, ... ,11. 

the k = 11 curve in figure 8 is given by 

Ao < 2(1 + e) + 2.683 6^ + ... (e < 1), (3.10) 

where d = 1 + e. Note that known superconformal theories populate the entire factorization 
hne/° so it is impossible to have a bound stronger than Aq < 2d. Our bound on dim($^$) is 
one of the few examples computed to date that approaches the provably best possible bound 
for some nontrivial range of d's. 

Eq. (3.10) can be directly tested in theories that admit a perturbative Banks-Zaks limit 
and contain a chiral operator with dimension near 1. As far as we are aware, there are 
no known examples of perturbative theories living above the factorization line. Here we 
have shown numerically that this can be understood purely from the constraints of crossing 
symmetry and unitarity. It would be very interesting to understand this fact analytically. 

It is amusing to speculate on the form of the bound as — oo. A simple and intriguing 
possibility is that the small-d behavior might extend to all d, so that the best possible bound 
Ao < 2d is realized. In other words, it might be the case that the anomalous dimension 

^"^Namely supersymmetric mean field theories, which satisfy the necessary requirements of unitarity and 
crossing symmetry, and exist for each d > 1. They occur in the infinite- iV limit of supersymmetric gauge 
theories. 
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Upper bound on dim($^$) 




Figure 7: An upper bound on the dimension of <I>1^<I>, where <^ is a chiral primary scalar of dimension 
d in an SCFT. The dashed hne is the factorization value A = 2d. Here we show k = 2, ... ,11. 

7$t$ = dim($'''<l>)— 2 dim($) is always non-positive. This possibility was investigated recently 
for theories with a weakly-coupled gravity dual in [56], with inconclusive results; effective 
field theories in AdSs allow for both positive and negative contributions to 7$t<i.. However, 
it's possible that additional constraints might be present in those theories which admit a 
consistent UV completion. 

Another possibility is that the bound converges above the factorization line, with a shape 
similar to the k = 11 curve in figure 7. In that case, one might wonder about the significance 
of the cusp near d = 1.4, which appears to be a common feature of each curve with k > A. A 
previous example of a dimension bound with a cusp is the 2D real scalar dimension bound, 
presented in [50] (building on the first 2D results of [49]). There, an actual theory, the 2D 
Ising model, exists very near the cusp, so that the bound is close to the best possible at that 
value of d. By analogy, one might speculate that an A/" = 1 SUSY 'minimal model' exists in 
the cusp in figure 7. 

3.4.1 Phenomenological Applications 

Our bound on dim($^^<l>) has implications for several models that use strong superconformal 
dynamics to tailor soft parameters in the MSSM. One example is the solution to the fi/Bfi 
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Upper bound on dim($^$) 
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Figure 8: A zoom in on the region of figure 7 near dim($) = 1. 



problem in gauge mediation proposed in [31, 32] and further developed in [33-36]. In this 
scenario, SUSY breaking is communicated to the visible sector via a chiral field X which 
develops a SUSY-breaking VEV {X) = FO"^ at some scale Air. In matching to the MSSM 
at Air,, the effective operators 



O 



X 



+ h.c. 



and 



^ X^ XHuHd 



M2 



(3.11) 



contribute to /i and -B/i, respectively. Here, M^, is the scale where these operators originate 
(typically the messenger scale). Many of the simplest gauge- mediated models generate both 



Ox and C^tx at one-loop at the messenger scale, so that naively cx ~ c^tx 



167r2 ' 



with A an 



0(1) coupling constant. However, this then leads to the problematic relation Bfi/ fi"^ ~ IGtt^, 
which precludes viable electroweak symmetry breaking. 

The solution proposed in [31, 32] is that X should participate in strong conformal 
dynamics over some range of scales Air < fx < Auv, with Auv < M*. If the anomalous 
dimension 7xtx = dim(X"''X) — 2 dim(X) is positive, then the operator Ox^x will be 
suppressed relative to Ox, and B^/fi"^ can be close to unity at the matching scale Air. 
In particular, to restore proper electroweak symmetry breaking, we should approximately 
have 



/AiRV-tx 

Uuv; 



167r2 



(3.12) 
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Using this relation, our upper bound on dim(X^X) in figure 7 translates into a lower bound 
on the running distance Ayy/Aj^, shown in figure 9. Note in particular that a small dim(X) 
requires a very large running distance, since our bound on 'jx^x approaches zero as dim(X) — )■ 
1. Consequently, viable models should at least have dim(X) > 1.3. Note that dim(X) can 
almost always be calculated using a-maximization in concrete examples, so a bound on the 
required running distance can be easily read from figure 9 for specific models. 

Running distance needed to solve fJ^/Bfi 




J — ^ ^ ^ , ^ ^ ^ , ^ ^ ^ , ^ ^ ^ , ^ ^ 

1 1.2 1.4 1.6 1.8 

Figure 9: An approximate lower bound on the running distance required for solving the l-i/Bfi 
problem with strong conformal dynamics, as a function of (i = dim(X). The middle curve 
corresponds to a loop factor suppression: Cxtxi^iR) = y^c^(^ir)) while the outer curves 
correspond to suppressions within factors of 2 and 5 of a loop factor. 

Our bound can also apply to models of conformal sequestering [19, 20, 22-27, 32] which 
contain chiral gauge singlets, where the idea is that a large dim(X'''X) can lead to suppression 
of fiavor-dependent soft-mass operators, 

[ d'e^X^x4<f>,. (3.13) 

Let us for example assume a gravity mediated scenario, where the cutoff scale is ~ Mpi 
and conformal running occurs between Mpj and an intermediate scale Aint ~ 10^^ GeV. 
Viable fiavor physics then roughly requires dim(X''"X) — 2 > 1 [26], and from figure 7 we see 
that such models should also have dim(X) > 1.35 or so.^^ Our bounds similarly constrain 

^^However, it's possible that one could avoid these constraints by having 'safe' flavor currents appear in 
the OPE (as discussed in [26]). 
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the possible suppression of these operators in superconformal flavor models [14-18, 21, 28- 
30], where the visible sector fields participate in the strong conformal dynamics. Once again, 
in all of these situations a comparison to our bounds can be checked in concrete examples 
using a-maximization. 

4 Bounds on OPE Coefficients 

In this section we will turn our attention away from bounding operator dimensions and 
instead explore some of the more basic bounds on OPE coefficients obtainable using these 
methods. We'll begin by reproducing (and strengthening) the upper bounds on scalar OPE 
coefficients for general CFTs previously presented in [51]. Then we'll focus on something 
qualitatively new — the possibility of placing lower bounds on OPE coefficients in theories 
that have a gap in the spectrum of operator dimensions. In fact, this happens naturally in 
supersymmetric theories for protected operators appearing in the $ x $ OPE, where a gap 
is forced by unitarity. We will then demonstrate that there are extremely constraining upper 
and lower bounds on the OPE coefficients of these operators when dim($) < 3/2. 

4.1 Scalar Operators in General Theories 

Let us begin by producing bounds on OPE coefficients of scalar operators Oq of dimension 
Aq appearing in the 0x0 OPE, where is a scalar operator of dimension d. As we saw 
in Eq. (2.39), by applying a linear functional a to the CFT crossing relations we can obtain 
an upper bound A^^^ < — a(Fo,o)- In figure 10 we show the best upper bounds on Xoq as a 
function of Aq that we have obtained so far, for d = 1.01, . . . , 1.66 with a spacing of 0.05. 
These bounds are obtained using k = 11, corresponding to a 6 6- dimensional search space. 
This plot strengthens bounds previously presented in [51]. 

Figure 10 clearly contains a lot of interesting structure. First, as (i — )■ 1, the curve becomes 
more and more sharply peaked around Aq — 2, with the height of the peak converging to 
the free value Aq = a/2.^^ On the other hand, as Aq — t- 1 all of the curves drop sharply 
to zero (first peaking at larger values of d), corresponding to the fact that a free operator 
cannot appear in the OPE. All of the bounds also increase in strength as Aq becomes large, 
possibly asymptoting to zero. Finally, as d increases at fixed Aq the bounds monotonically 
decrease in strength. Note that in the present study we have found the region d > 1.66 
to be numerically more difficult (though very weak bounds appear to exist at least up to 
d ~ 1.86), and we postpone a full investigation of this region to future work. 

Let us take a moment to understand a way in which our method fails to fully pick out the 
spectrum of free theories as d ^ 1. While our upper bound becomes nicely peaked around 
the free value in this limit, our algorithm cannot easily distinguish between a single Aq — 2 
operator with Aq — \/2, and a broader spectrum of operators, each having Aq somewhat 

^^Note that the free OPE coefhcient is \/2 rather than 1 because we have required the 0^ operator to have a 
canonically normahzed two-point function, rather than the normahzation inherited from Wick contractions. 
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close to 2 and Aq < "\/2- The issue is that both of these scenarios can lead to very similar 
conformal block contributions to the 4-point functions that we are studying. On the other 
hand, if we knew that there was only a single operator appearing in the OPE up to a certain 
dimension, this ambiguity could not occur and we would be able to also place lower bounds 
on its OPE coefficient. In the next subsection we will study this possibility in more detail, 
focusing on protected operators appearing in the $ x $ OPE in SCFTs. 



Upper bounds on scalar OPE coefficients, d = 1.01,. ..,1.66 




Figure 10: Upper bounds on the OPE coefficient of a scalar operator Oq € (p x (p (not necessarily 
of lowest dimension). Each curve is for a different value d = 1.01, . . . , 1.66, with a spacing of 0.05 
and d = 1.01 corresponding to the lowest curve. Here we have taken k = 11. 



4.2 Protected Operators in Superconformal Theories 

As we reviewed in section 2.3, if $ is a chiral superconformal primary of dimension d in an 
Af = 1 SOFT, the $ X OPE contains superconformal primaries of dimension A > i + 2 
and their descendants. On the other hand, the $ x $ OPE can contain a chiral $^ operator 
of dimension 2d, superconformal descendants QO^ of protected operators having dimension 

— 2 

2d + i, and superconformal descendants Q O oi unprotected operators with a dimension 
satisfying A > \2d - 3\ + 3 + i. 

Notice that, as long as d < 3/2, there is necessarily a gap between the dimensions of 
the protected operators appearing in the $ x $ OPE and the dimensions of the unprotected 
operators. This gap is a consequence of the unitarity constraints on operator dimensions 
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in SCFTs. Because of this gap, no other operators appearing in the OPE can give similar 
conformal block contributions to the four-point function ($$^$$^), and we can attempt to 
derive /ot^er bounds on the OPE coefficients A<j>2 and Ag^^, in addition to upper bounds. 

The logic used to obtain a lower bound requires only a slight modification to the procedure 
described in section 2. Since one could in principle attempt to obtain a lower bound in any 
theory with a dimension gap, let us first describe the logic for the simplest case of the real 
scalar crossing relation in general CFTs. To obtain a lower bound on an OPE coefficient 
-^Oo' '"^^ again consider applying a linear functional to the real scalar crossing relation, 
as in Eq. (2.36). However, instead of imposing the constraints (2.37) and (2.38), we can 
alternatively require 



a(^Ao,£o) = 1' and 
(^{F^Ai) < 0, for all other operators in the spectrum. 



(4.1) 
(4.2) 



which leads to the lower bound 

-a(^o,o 



\2 



E 



pos. X neg. > -a(Fo,o). 



(4.3) 



Note that (4.1) and (4.2) are only compatible with each other if we know that there is a gap 
between Aq and the A's for all other operators in the spectrum. 

Generalizing to the superconformal crossing relation of Eq. (2.31), if we isolate a protected 
operator Oq of spin Iq and require 




< 0, 



for all other operators in $ x $, and 



< 0, for all (non-unit) operators in $ x <l>^. 



(4.4) 



(4.5) 



(4.6) 



we obtain the lower bound 



a: 



do 



> 




(4.7) 



Meanwhile, reversing the inequalities in (4.5) and (4.6) leads to an upper bound on A^,^, 
following our usual logic. 

In figure 11 we show the resulting upper and lower bounds on A$2, where we have taken 
k = 2, . . . , 11 in the numerical optimization. We can see that the strongest bounds are 
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extremely constraining when d = dim($) is even somewhat close to 1, forcing A$2 to live 
very close to the free value A<j,2 = a/2. In particular, these results imply that it should not 
be possible to construct a weakly-coupled (Banks-Zaks) SCFT where both d and A<i,2 are 
modified at the one-loop level. Indeed, in all constructible examples A$2 receives its leading 
correction at second order in perturbation theory. On the other hand, we see that the lower 
bound disappears before d = 3/2, as expected, while the upper bound persists. 

As (i — )■ 2, we may also compare the upper bound to the OPE coefficients of composite 
operators in theories containing free chiral superfields. In the simplest case, we can consider 
a single free field Q and then identify $ = -^Q"^- In this case the operator $^ = -^Q^ is 
canonically normalized, so the OPE is 

^-:Q') x(^^Q')-^(4^Q']+---^ (4-8) 



V2 J \V2 

and we have A<j,2 = -^/G, which is consistent with the bound. More generally, considering the 
dimension-n operator $ = "TtQ" leads to an OPE coefficient of A$2 = -^^^^^l — , which the 
bound must respect at even higher integer values of d. 

Another simple generalization is to consider meson operators M = —^Q^Qi built out of 

N free quarks Q^. In this case Wick contractions give a two-point function ((M^)(M^)"'') ~ 
2 + so the OPE in terms of canonically normalized operators is given by 



M ^ M J9. + — \ — + (4.9) 




Thus, we can read off an OPE coefficient of A$2 = ^y2 + A/N, which is consistent with our 
bound for all values of A^. It is interesting to see that while OPE coefficients of composite 
operators with c? ~ 2 know about the underlying constituents of the operator, as (i — )• 1 
the OPE coefficient necessarily loses memory of where the operator came from. Indeed, free 
operators have no hair! 

In figure 12 we extend these upper and lower bounds to OPE coefficients of the other 
protected operators QOe appearing in the $ x $ OPE. Here we give the results for i = 
2, . . . , 10 and have taken A; = 11 in the numerical optimization (though similar bounds also 
exist at larger values of i). All of the bounds continuously interpolate to the free values as 
(i — >■ 1, given by Xqo^ = ^j^^TTi- Notice that all lower bounds vanish before c? = 3/2, as 
they should. ^'^ 

Taken together, the upper and lower bounds on Ag^^ are extremely strong, almost 
determining this coefficient when ci < 1.4. One can view this singling out of an essentially 
unique OPE coefficient as a remarkable success of the 4D conformal bootstrap program! It 
is worth comparing the bounds to the known values of Ag^^ in supersymmetric mean field 



^•^Once they are computed, one can include information about these lower bounds in semidefinite programs 
for other quantities, like e.g. upper bounds on OPE coefficients of operators in the $ x <I>'f OPE. We found 
that this procedure does not significantly improve the results in practice. 



33 



Upper and lower bounds on A<i,2 

A<J)2 




Figure 11: Upper and lower bounds on the OPE coefficient of <1>^ in <^ x <I>, as a function of 
d = dim($). The dashed line indicates the free value A$2 = \/2. The points shown at d = 2 

indicate the sequence of values A$2 = \J2 + realized for composite operators in free theories. 
We give the bounds for A; = 2, ... , 11. 

theories (MFTs), which occur in the planar limit of large- gauge theories. There, the role 
of QOi is played by the 'twist-2(i' double-trace operators 

Of^ = <l>cf^i---cf'^^<l>- traces, (4.10) 

with even spin i. Their (squared) OPE coefficients in $ x $ are given by [68] 

2 _ 2T\d + i)T{2d + i -1) 
of^ ~ V\d)V{l + l)V{2d + 2t-iy 

and these values of A (2) are shown as dashed lines in figure 12, for £ = 2, 4, . . . , 10. They are 
fully consistent with both our upper and lower bounds on Ag^^. Note that the MFT value 
of A (2) is equal to the free value -\/2, so it is consistent with our bounds in figure 11. 

The striking agreement between our bounds and the mean field theory values of OPE 
coefficients at small d has interesting implications for SCFTs with weakly-coupled AdSs 
duals. In such theories, corrections to OPE coefficients away from their MFT values can be 
computed in perturbation theory using Witten diagrams. Our bounds imply that corrections 
to A^ must vanish to very high order in [d — 1), particularly at large £. If any corrections 

Of 
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(4.11) 



0.001 



Upper and lower bounds on Ag^^, £ = 2, 4, . . . , 10 




d 



Figure 12: Upper and lower bounds on the OPE coefficients of protected operators QOi appearing 
in $ X along with their mean field theory values Eq. (4.11) (dashed lines), for ^ = 2, 4, . . . , 10. 
Each curve goes continuously to the free value V^ (2^)i'i/2 at d = 1. All lower bounds vanish at 

d = 3/2, since the gap in dimensions between QOn and non- protected operators disappears at that 
point. Here we have taken k = 11. 



were nonzero at finite values of (d — 1), then we would obtain sharp bounds on bulk coupling 
constants. We defer further exploration of these interesting constraints to future work. 



5 Bounds on Central Charges 

In this section we explore bounds on the OPE coefficient appearing in front of the stress 
tensor T^^ , which is a conserved spin-2 operator of dimension 4 that must be present in any 
CFT. Since this OPE coefficient is fixed by a Ward identity in terms of the central charge c 
of the theory (defined as the coefficient appearing in the two-point function (T^^'^T^^) oc c), 
we will ultimately be deriving bounds on c. Previously, lower bounds on the central charge 
in both general CFTs and SCFTs were explored in [52, 53, 55]. The main new results of 
this section will be to extend these analyses to situations with global symmetries, where we 
will show that there are bounds on the central charge that scale with the size of the global 
symmetry representation. 
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5.1 General Theories 



Let us begin by establishing some notation. The stress tensor is typically normalized as 

(r-n.)r^'(o)) ^ (5.1) 



where I^"'{x) = rj^'^ — and c is the central charge appearing in the trace anomaly, 

(T^) = j^(Weyl)^ — (Euler), when the theory is placed on a curved background. In 
this normalization a free scalar has Cfrcc = and a free Weyl fermion has Cfreefermion = ^• 

The stress tensor is the local current generating the dilatation charge, where in ra- 
dial quantization D = — J dQ x^XyT^'-^ (the integral is over a three-sphere surrounding 
the origin). Requiring the action -D0(O) = (i0(O) then fixes the OPE to have the form 
T'^'^(x)0(O) ~ — 3^2^,6 {x^x'^ — ■^Tj^'^x'^^ 0(0) "!"•••) which leads to the stress tensor conformal 
block contribution 

a^?2a^34(##) ~ ^^74,2 (general CFTs). (5.2) 

Generalizing to the situation where 0* transforms under an SO(A^) or SU(A^) global sym- 
metry, the stress tensor appears as an operator in the sum rules given in Eqs. (2.14) 
and (2.23), again with OPE coefficient = Note that a free real scalar transforming 
as an SO(A^) fundamental or a complex scalar transforming as an SU(A^/2) fundamental 
gives a contribution of A^Cfree to the central charge. 

To begin, in figure 13 we show the bounds on c obtained by applying our semidefinite 
programming algorithm to the case of a single real scalar 0, where we show curves for 
/c = 2, . . . , 11 in the numerical optimization. We see that for k > Q, the bounds smoothly 
approach the free value Cfree as d — )■ 1. This is consistent with and improves upon the 
bounds on c previously presented in [52, 53]. Note that here we are only assuming that the 
dimensions of operators appearing in the 0x0 OPE satisfy the unitarity bound — one could 
also assume that is the lowest dimension scalar in the theory to obtain somewhat stronger 
bounds at larger values of d as was done in [53]. However, here we make only the minimal 
assumption to allow for a more straightforward comparison to our other bounds. 

In figure 14 we show bounds on c in the presence of SO(iV) or SU(iV/2) global symmetries 
for = 2, . . . , 14. Here have taken k = 11. We see that the bounds also smoothly approach 
the free values as d — )■ 1, scaling linearly with N. This greatly improves upon the bounds 
derived in [52, 53] (and given in figure 13) for theories with global symmetries. The reason 
for the improvement is that here we have incorporated the constraints of crossing symmetry 
for all operators in the 0* multiplet; without doing this the bounds of [52, 53] could not 
differentiate between the stress tensor and other spin 2 operators (e.g., the SO(iV) symmetric 
tensor (p^'^d^d'^cp^^) that have an 0(1) OPE coefficient in the d — )■ 1 limit. 

It is interesting to understand the implications of the bound of figure 14 for the AdS / CFT 
correspondence. For theories with an AdSs dual description, the bulk Planck scale is 
proportional to c, the bulk gauge group is identified with the SO(A^) or SU(A^/2) global 
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Lower bound on c for a real scalar 

c 

0.016 r 




Figure 13: A lower bound on the central charge of a theory containing a scalar (j) of dimension 
d. The dashed line indicates the value Cfree = 1/120, corresponding to the central charge of a free 
scalar. Here we show bounds for the values k = 2, ... ,11. 

symmetry, and d is related to the masses of bulk fields. Our bound then says that theories 
with sufficiently light bulk excitations cannot have a gravitational scale that is arbitrarily 
small. Moreover, if those fields transform as fundamentals under the bulk SO(A^) or SU(A^/2) 
gauge group (and correspond to operators with ~ 1), then the Planck scale must scale at 
least linearly with N. 

It would be fascinating to identify CFTs that live close to these bounds, particularly 
in the large N limit. Unfortunately, in gauge theories believed to flow to conformal fixed 
points that also posses an SO(A^) or SU(A^/2) global symmetry, the central charge typically 
scales as A^^ , at least near (i ~ 1 . The reason is that conformality forces the size of the global 
symmetry to scale proportionally to the size of the gauge group, and gauge degrees of freedom 
live in adjoint representations of the gauge group which have 0{N^) components. We will 
see examples of this in the next subsection, where we extend the bounds to superconformal 
theories in which c is explicitly calculable. 
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Lower bounds on c for SO(iV) or SU(iV/2), iV = 2, . . . , 14 

c 




Figure 14: A lower bound on the central charge of a theory containing a scalar (pi of dimension d 
transforming as a fundamental of an SO(A^) or SU(A^/2) global symmetry, for = 2, . . . , 14. In 
this plot Cfrcc = 1/120, corresponding to the central charge of a free scalar. Here we have taken 
k = 11. 

5.2 Superconformal Theories 

In A/" = 1 SCFTs, the stress tensor is a superconformal descendant of the spin-1 U{1)r 
current, T ~ {QQJr)i+i, as in Eq. (2.25). Applying Eq. (2.26) to (5.2), we see that has 
an OPE coefficient oi \\ = appearing as an S'^ operator in the superconformal sum 
rules of Eqs. (2.31) and (2.33). Since a free chiral superfield contains both a complex scalar 
and a Weyl fermion, it gives a contribution of Cchirai = 2x^1^ + ^ = ^. 

In figure 15 we show the results of our semidefinite programming algorithm for obtaining 
bounds on the central charge of any theory containing a chiral scalar $. We give the results 
for k = 2, . . . , 11, where all of the curves for k > 3 drop sharply very close to (i ~ 1 and 
go just below the free value. The k = 11 curve significantly improves upon SCFT central 
charge bounds previously obtained in [52, 55]. Note that the sharpness of the drop (reaching 
within 1% of the free chiral value closer than d ~ 1.0000002) is strong evidence that the free 
theory is an isolated solution to the crossing relations. This is intuitive from the perspective 
of constructing perturbations of the free theory — all such perturbations leading to an 
interacting SCFT require additional matter, which increases the central charge. In order to 
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demonstrate that the bound does in fact approach the free value, in figure 16 we also show 
the bound for k = 11 where {d — 1) has been placed on a logarithmic scale. 



Lower bound on c for a chiral scalar 

c 




Figure 15: A lower bound on the central charge of any SOFT containing a chiral scalar $ of 
dimension d. The dashed line is at Cchirai = 1/24, corresponding to the central charge of a free 
chiral superfield {d = 1). Despite appearances at this zoom level, all the curves above drop sharply 
near d = 1 and interpolate smoothly to the free value. In this plot we have taken k = 2, ... ,11. 

We extend these bounds to the situation where $* transforms as a fundamental under 
an SU(A^) global symmetry in figure 17, where we have taken k = 10 and show curves for 

= 2, . . . , 14. All the curves interpolate to the free values A^Cchirai as — )• 1, in all cases 
with a very sharp drop in the bound close to 1. Again we see that the bounds scale linearly 
with A^, and moreover the linear behavior extends out to larger values of d compared to the 
non-supersymmetric bounds of figure 14. 

Let us now take a moment to compare these bounds to some concrete SCFTs. The 
reason that such a comparison is possible is that both d and c are calculable in terms of the 
U{1)r symmetry — d is calculable because the dimensions of chiral superconformal primary 
operators are related to their R charge a.s d = |-R, and c is calculable via 't Hooft anomaly 
matching using the relation c = ^(9Tri?^ — 5Tri?) [69, 70]. The U{1)r symmetry can then 
often be determined using symmetry arguments, or more generally using a-maximization [71]. 

One of the simplest = 1 SCFTs is supersymmetric QCD with gauge group SU(A''c) 
and Nf flavors of quarks Q, Q in the conformal window ^Nc < Nf < 3Nc [57]. In this case 
the gauge-invariant mesons M = QQ have dM = 3(1 — Nc/Nf), while the central charge 
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Lower bound on c for a chiral scalar 



c 




Figure 16: The k = 11 curve of figure 15, where {d — 1) has been placed on a logarithmic scale. 
The bound smoothly approaches the free value Cchirai = 1/24 very close to d = 1. 



is evaluated as c = j^ilN"^ — dN^/Nj — 2). The mesons are bi-fundamentals under the 
SU(A''/) X SU(A''/) symmetry group, so our bounds will apply by considering either of these 
groups. 

However, we immediately see that the central charge in SQCD grows like 0(A^^), so 
theories at large values of A'^^ ~ trivially satisfy the bounds. On the other hand, all of 
the small theories still have a central charge larger than 1 = 24cchirah so the bound is 
also easily satisfied for these theories. Part of the problem is that we have only included 
a subgroup of the full SU(A^/) x SU(A^/) global symmetry when deriving our bounds. In a 
future publication [72] we hope to extend the bounds to bi-fundamentals transforming under 
an SU(A^) x SU(A^) symmetry group, in order to make closer contact with the values realized 
in SQCD and similar theories. 



6 Bounds on Current Two-point Functions 



6.1 General Theories 



Now let us turn to placing bounds on another set of fundamental OPE coefficients, namely 
those appearing in front of spin-1 conserved global symmetry currents. In the OPE between 
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Lower bounds on c for a SUSY SU(A^) chiral scalar, = 2, . . . , 14 




d 



Figure 17: A lower bound on the central charge of any SOFT containing a chiral scalar <I>j of 
dimension d transforming as a fundamental of an SU(A^) global symmetry, for = 2, . . . , 14. Here 
Cchirai = 1/24 denotes the contribution to c from a free chiral superfield. Despite appearances at 
this zoom level, all the curves drop sharply very close to d = 1 and interpolate continuously to the 
free values. In this plot we have taken k = 10. 

SO(A^) or SU(A^) fundamentals, we should be careful to distinguish between the SO(A^) or 
SU(A^) symmetry currents living in the adjoint representation and singlet currents associated 
to some other global symmetry that we are not considering explicitly. 



6.1.1 Adjoint Currents 

Let us begin by focusing on the case of adjoint currents. Consider a CFT with some 
global symmetry, containing a scalar field (pi transforming in some representation of this 
symmetry. We will denote by T^^ the generators in this representation. The associated 
conserved currents transform as global symmetry adjoints. Ward identities completely fix 
the three-point functions with one current insertion: 



i(Xi)0j(x2)J;f(x3)) 



27r2 



2-2d 

^2 ^2 ' 
•^13-^23 



where Z,, 



X 



13 



•^12 



(6.1) 
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With the above normahzations, the two-point function (J'^J^) contains undetermined coef- 
ficients that roughly measure the amount of stuff charged under the global symmetry: 



(J^-(xOJ^^(x,)) = ^^^^. (6.2) 

Let us write r^'^ = kTt{T^T^), where k can be viewed as a symmetry current 'central 
charge.' As we did for the energy momentum tensor, we can rescale to have a canonically 
normalized two-point function and absorb k into the OPE coefficient Aj associated with the 
current. In the end, the contribution of an adjoint current to a four-point function of 0j's 
can be written 

xiixli{M,<l>k<l>i) ~ -^Tr(T^T^)-i7;X 93,1- (6.3) 

In order to proceed further we need to specify the global symmetry group. For instance, for 
SO(A^) and (pi in the vector representation, one can show that 

TT{T^T^rX^T,^i = li^uS.k - S^kSji), (6.4) 

and consequently, comparing to Eq. (2.13), we have Aj = -j^. Similarly, for SU(A^) and 0, 
in the fundamental representation, we have 

Tr(T^T^)-i(T^)5(T^)f = StS^-^S^S^, (6.5) 

so that Aj = g^. These relations hold for currents appearing in OPEs in general CFTs; 
we will discuss the generalization to A/" = 1 superconformal theories below. However, first 
we will consider the situation of singlet currents appearing in the OPE, namely currents 
corresponding to a global symmetry that is different from the SO(A^) or SU(A^) that we are 
studying. 



6.1.2 Singlet Currents 

As mentioned above, the SO(A^) or SU(A^) global symmetry current is not the only conserved 
spin-1 operator of dimension 3 that can contribute to the four-point function; additional 
currents, possibly transforming in different representations, may also exist. Clearly the 
presence of an additional conserved current implies the existence of a global symmetry beyond 
the one exploited to write the crossing symmetry constraints. The OPE coefficient associated 
to this operator not only contains the two-point function normalization, but also parametrizes 
our ignorance about the nature of the additional global symmetry. Indeed, when the global 
symmetry is not specified the three-point function coefficient could in principle be arbitrary. 

In the case of fundamentals transforming under an SO(A^) global symmetry, spin-1 opera- 
tors appearing in the OPE can only transform in the adjoint (antisymmetric) representation, 
corresponding to the SO(A^) current itself. In the case of SU(A^) fundamentals, along with 
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the adjoint current we also have the possibihty of SU(A^) singlet currents. For example, we 
can think about a CFT with a global symmetry SU(A^) x Q. If we consider scalar operators 
transforming in some representation of Q with generators T^, then the ^-current is a singlet 
with respect to SU(A^), and its contribution to the four-point function will be 

^Ti^T'^T^y'^TifTiZ 93,1 = ><%j6img3,i, (6.6) 

where Kg is the two-point function of the ^-current. Until we additionally specify the Q 
symmetry group and charges, this parameter is arbitrary. However, we can collectively define, 
by analogy with the adjoint current, an effective current two-point function normalization 
Kes = l/6Aj. We will place bounds on Kes when we give our results below. 



6.1.3 PVee Theory and Numerical Results 

To clarify the above discussion, let us analyze in detail the theory of free complex scalars, 
using only information about the SU(A^) global symmetry, which is contained in the larger 
S0(2A^) symmetry of the theory. The OPE 0j x 0-'"'' contains an adjoint current and a singlet 
current, both conserved: 

Jt ~ Js ~ (6.7) 

The conformal block decomposition of the scalar four-point function directly gives us the 
values of the singlet and adjoint OPE coefficients for spin-£ currents: 

AL = |^. -^ = 5. (6.8) 

As = ^^, -efr = y. (6.9) 

The first point that we notice is the different scaling of the two above quantities with the 
size of the symmetry group. While the adjoint current two-point function normalization 
is independent of A^, the singlet one grows with the dimension of the representation. We 
therefore expect lower bounds on k^s to scale with A^, similarly to the way that the central 
charge bounds did in the previous section. 

Let us now discuss the same theory, using the whole SO(A^) global symmetry. This 
time only the adjoint current contributes to the four-point function, and its OPE coefficient 
(along with the other spin-£ adjoint operators) can be determined from the conformal block 
decomposition (see for instance [54]) as 

(£1)2 \ 



^'^In addition, the OPE x could contain conserved spin-1 operators transforming in the antisymmetric 
representation of SU(A^). However, such currents (along with their complex conjugates) would generate 
charges which enhance SU(7V) to a larger group SU(A^) — >■ S0(2A^). Thus, such theories necessarily fall 
under the class of CFTs with a global S0(2Af) symmetry, which we consider separately. 
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Now that we have an intuition for the free values of k and Kcs, we are ready to present 
numerical bounds in several classes of theories. In figure 18, we show a lower bound on the 
two-point function coefficient k, for a CFT with an SO(A^) global symmetry for = 2, . . . , 14. 
As expected, when ti — )■ 1, all of the bounds drop sharply to the free SO(A^) value k = 1/6. 
The bounds get stronger as A*" increases, while as d varies away from 1, they first become 
stronger and then weaken. 

Lower bounds on k, for SO(A^) adjoint currents, A^ = 2, . . . , 14 



0.25 



0.15 




0.05 - 



d 



Figure 18: A lower bound on the two-point function coefficient {Jf^J^) oc kTi{T^T^) of the SO(A^) 
adjoint current appearing m. (p x (j)^ where (j) transforms in the fundamental of an SO(A^) global 
symmetry group, for A^ = 2, . . . , 14. All curves smoothly approach the free SO(A^) value k, = 1/6. 
Here we have taken k = 11. 

As a second example, in figure 19 we consider the case of an SU(A^) global symmetry 
and present lower bounds on for a singlet current. Our expectation that the constraints 
scale almost linearly with A^ (when d is close to 1) is confirmed. Thus, this quantity serves 
as a rough measure of the number of degrees of freedom in the theory transforming under 
the symmetry, at least near d = 1. One the other hand, the linear scaling disappears as d 
increases. 



6.2 Superconformal Theories 

Let us generalize the above bounds to theories with M = 1 supersymmetry, where currents 
are descendants of scalar superconformal primaries of dimension 2. Consider four-point 
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Lower bounds on k^s for SU(A^) singlet currents, = 2, . . . , 14 




d 



Figure 19: A lower bound on the effective two-point function coefficient n^g = l/6Aj of SU(A^) 
singlet currents appearing in x 0^^, where (^j transforms in the fundamental of an SU(A^) global 
symmetry group, for N = 2, ... ,14. All curves interpolate continuously to the free values A^Kfree 
where K{ree = 1/3, and in this plot we have taken A; = 11. 



functions of chiral and anti-chiral operators transforming under an SU(A^) 

global symmetry. SU(A^) adjoint currents give a superconformal block contribution 



0) 



while SU(A^) singlet currents give an effective superconformal block contribution 

1 



x?^x^f($,$5t$^,$It) ^ —616102,0. 



(6.11) 



(6.12) 



In figure 20, we show bounds on n for adjoint currents appearing in $j x for SCFTs 
with an SU(A^) global symmetry and = 2, . . . , 14. These bounds again increase strongly 
with A^, growing as a roughly affine function. For d < 1.5, k must be substantially higher 
than its free value, with the bound dropping sharply to the contribution of a free chiral 
superfield Kchirai = 1 near d = 1. Consequently, the free theory appears to be isolated in 
the space of SCFTs with an SU(A^) flavor symmetry. This accords with our intuition from 
theories with a Lagrangian description. To couple a free SU(A^) fundamental to a nontrivial 
interacting sector (and thus raise its dimension away from d = 1), we need additional matter 
which must itself transform under SU(A^). 
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Lower bounds on k, for SUSY SU(A^) adjoint currents, = 2, . . . , 14 
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Figure 20: A lower bound on the two-point function coefficient (J^J^) oc kTt{T^T^) of an SU(A^) 
adjoint current appearing in $j x , where is a chiral scalar transforming in the fundamental 
of an SU(A'") global symmetry group in an SCFT, for = 2, . . . , 14. Despite appearances at this 
zoom level, all the curves above drop sharply near d = 1 and interpolate continuously to the free 
value Kchirai = 1- Here we have taken k = 10. 

In figure 21, we also show a lower bound on k^s for singlet currents appearing in $j x 
Once again, we see that these bounds increase with N, scaling roughly linearly for small d. 
As in the adjoint case above, the bounds drop very sharply to their free values A^^chirai near 
d = 1, while the scaling disappears as d increases. 

6.2.1 Comparison to SQCD 

As with central charges, our bounds on current two-point functions can be checked explicitly 
in a given superconformal theory. For example, in SUSY QCD, S\J{Nf)L and SU {Nf)^ flavor 
currents appear in the OPE of a chiral meson and its conjugate 



Here, i,j are indices for SU {Nf)^ and?, j are indices for S\J {Nf)^. We have not yet generated 
bounds that exploit the full SU(A^/)i x SU(A^/)r symmetry group of SQCD. However, we 
can compare to our SU(A^) bounds by 'forgetting' one of the flavor groups, say S'U{Nf)fi, 
and examining the theory from the point of view of SU {Nf)L alone. Specifically, we shall set 



M~ X M/ 



(6.13) 
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Lower bounds on KcS for SUSY SU(A^) singlet currents, = 2, . . . , 14 




d 



Figure 21: A lower bound on the effective two-point function coefficient k^s = 1/Aj of SU(A^) 
singlet currents appearing in $j x where <I>i is a chiral scalar transforming in the fundamental 
of an SU(A^) global symmetry group, for = 2, . . . , 14. Despite appearances at this zoom level, all 
the curves above drop sharply near d = 1 and interpolate continuously to the free values A^Kchiral 
where Kchirai = 1- Here we have taken k = W. 



t = J = 1, so that the right-flavor currents are then singlet scalars in 
the left-flavor currents are adjoints. 



X M/, while 



The current two-point functions for Jjj and in SQCD both scale like Nf (or Nc). 
However, only our SU(A^)-singlet bounds scale with A^, and thus have a chance of approaching 
the values for SQCD. Consequently, we will focus on the contribution of to the conformal 
block expansion of meson four-point functions. This reads 



xlixliiMl'Mj^Ml^Mi^) 



(6.14) 

.AB\~l 



where ^2,0 is the superconformal block for a conserved current multiplet and tab = [t 
is the inverse two-point function coefficient for J^. In superconformal theories, t^^ can be 
computed simply in terms of 't Hooft anomalies using t^^ = — 3Tr(i?T"^T^ 
becomes 



For J^, this 



\JrJr) 



oc r 



AB 



^Nl^AB 



2N 



f 



where the SU(A^/) generators are normalized according to Tx(T^T^^ 



1§AB 



1 



3N^ 



6l6l 



N. 



(6.15) 
Thus, we have 
(6.16) 
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Lower bounds on KcS for SUSY SU(A^) singlet currents and comparison to SQCD 
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Figure 22: A lower bound on the effective two-point function coefficient Kcfj 



1/A2 of SV{N) 



singlet currents appearing in $j x where <I>i is a chiral scalar transforming in the fundamental 
of an SU(A'^) global symmetry group, for = 2, . . . , 14. Here we have taken k = 10. We have also 
plotted points corresponding to SQCD theories with various values of Nj and Nc- The lines below 
each point indicate the distance to the corresponding bound. Many SQCD theories lie within an 
0(1) factor from our bounds. 

In figure 22 we compare this value of KcS for several SQCD theories to our singlet current 
bounds from figure 21. For many values of Nf and Nc, our bound comes within an 0(1) 
factor of the SQCD value, with the smallest separation at small dimensions c? ~ 1. We 
expect our bound to become stronger with the added information of Sl]{Nf)ii symmetry, 
perhaps resulting in a hybrid of figures 21 and 20. It will be interesting to compare SQCD 
to these new bounds, and understand more about the structure of four-point functions in 
this important theory. 



7 Conclusions 



Let us briefly summarize our main results. In this work we explored bounds on operator 
dimensions and OPE coefficients in 4D CFTs and = 1 SCFTs, building on the previous 
studies performed in [49-55]. These bounds can be viewed as the initial stages of a concrete 
implementation of a 4D conformal bootstrap program. Here we focused on bounds in the 
presence of SO(A^) and SU(A^/2) global symmetries, which had previously shown themselves 
to be more difficult (but not impossible [55]) to obtain using algorithms based on linear 
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programming methods. In order to push the program further, we presented a new algorithm 
based on semidefinite programming, which utihzed the fact that derivatives of conformal 
blocks can be arbitrarily well approximated by positive functions times polynomials in the 
operator dimensions. This new algorithm enabled us to show that there are completely gen- 
eral bounds on CFTs and SCFTs in the presence of global symmetries that are significantly 
stronger than were previously known to exist. 

In particular, we greatly strengthened bounds on dimensions of singlet operators ap- 
pearing in the OPE between fundamentals transforming under SO(A^) or SU(A^/2) global 
symmetries. Bounds on dimensions of singlet operators in the presence of SO (4) or SU(2) 
global symmetries are relevant for models of conformal technicolor, and our bounds place 
severe constraints on these models, particularly when one does not assume any special flavor 
structure in four-fermion operators. We refer readers to [65] for further discussion of these 
constraints. In fact, in the present work we saw that bounds on singlet operators were in 
general identical between SO(A^) and SU(A^/2) global symmetries. We have so far not been 
able to construct a rigorous proof of this equivalence, so it would be good to gain a better 
understanding of it in future work. 

We also obtained similar bounds on operator dimensions in A/" = 1 SCFTs, where we 
showed that there are bounds on the lowest- dimension scalar appearing in the $ x OPE 
that appear to asymptote to the line A = 2d near (i ~ 1. This result is particularly interesting 
in light of the discussion of [56] on positive anomalous dimensions of these operator in SCFTs 
— our results demonstrate that this should not be possible when one is sufficiently close to 
the free limit. 

In this work we also initiated an exploration of both upper and lower bounds on OPE 
coefficients of protected operators appearing in the $ x $ OPE in SCFTs. In this case, 
lower bounds are possible due to the fact that there is a gap in the dimensions of operators 
appearing in this OPE that is required by unitarity. Because one can obtain bounds in both 
directions, we are able to see that the possible behavior is very tightly constrained even when 
one is only somewhat close to the free limit. We expect that similar lower bounds should be 
possible in any situation (including non-supersymmetric theories) where one assumes that 
there is a dimension gap such that only a single operator can contribute to the conformal 
block decomposition up to a certain dimension. 

We also explored bounds on central charges and current two-point function coefficients 
in the presence of operators transforming as fundamentals under SO(A^) or SU(7V/2) global 
symmetries, finding bounds that scale linearly with when the operator dimension is close 
to 1. An exception is the case when the current is the adjoint current corresponding to the 
SO(A^) or SU(A^/2) symmetry itself, in which case the bounds approach a value independent 
of in the free limit. In superconformal theories, these bounds can be compared to concrete 
theories where the central charge c and current two-point functions k are calculable using 't 
Hooft anomaly matching. While the central charge bounds are still relatively far from their 
realized values, we showed concretely that our bounds on k are an 0(1) amount away from 
the values realized in supersymmetric QCD in the conformal window. 

A clear future direction is to generalize these bounds on A/" = 1 SCFTs to situations 
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with bi-fundamentals transforming under SU(A^) x SU(A^) global symmetries (or adjoints 
transforming under SU(7V) global symmetries). Then one would hope to see bounds on the 
central charge that scale like ~ A^^, as well as significantly stronger bounds on current two- 
point functions. It will be fascinating to see how these bounds compare to concrete Af = 1 
theories such as supersymmetric QCD in the conformal window, particularly if one can find 
theories that nearly saturate the bounds. One could also input all known information about 
these theories and attempt to find even stronger constraints on the dimensions of unprotected 
operators. We plan to explore these bounds in a future publication [72]. 

Another interesting direction would be to apply these methods to four-point functions of 
operators with spin, such as symmetry currents or the stress tensor. To do this, one needs 
a tractable way of working with higher-spin conformal blocks. Recently, some progress has 
been made in this direction [75, 76], though more work may be needed in order to make a 
completely general analysis possible. However, if this program could be carried out, one could 
for example start to study whether crossing symmetry of stress tensor four-point functions 
is connected to the bounds on a/c obtained in [77]. In fact, it may be more immediately 
tractable to begin such explorations for theories with M = 2 supersymmetry, where the 
stress tensor is contained in a scalar multiplet. Similarly, one can explore crossing symmetry 
of current four-point functions in A/" = 1 theories, where progress at deriving the relevant 
superconformal blocks was made recently in [61]. We believe that these directions may be 
worth pursuing in future work. 

One would additionally like to generalize these bounds to 3D CFTs, where progress has 
been recently made at understanding the properties of 3D conformal blocks [78]. While 
closed-form expressions are not yet known, it is likely that recursion relations similar to 
what we used in the present study could make a numerical exploration tractable. It would 
be particularly interesting to see if one could learn more about the 3D Ising model using 
these methods, or whether bounds could be placed on the behavior of real-world condensed 
matter systems. It is also interesting to explore these bounds in 2D (expanding on the 
preliminary studies of [49, 50]) or in 6D, where perhaps progress can be made at unraveling 
the structure of the mysterious 6D (2, 0) SCFTs. 

Of course, it would be nice to have a better analytical understanding of the structure of 
the optimal bounds. While such an understanding has eluded us so far, it is possible that a 
new approach (such as studying the Mellin representation as in [79-81]) could shed light on 
the origin of these bounds. Less ambitiously, it would be good to study whether expansions 
of the crossing relation around other points in {z^ z) space may provide a more efficient way 
to find an optimal linear functional. A related question is to understand whether any of the 
multiple crossing relations that we have used in cases of global symmetries are redundant or 
unnecessary for obtaining an optimal bound. We leave such questions to future work. 

Finally, we hope that progress can be made at understanding where these bounds fit 
in the context of the AdS/CFT correspondence [38-40]. Bounds on the central charge 

^^An alternate approach to learning about these dimensions is to look for hidden structure such as 
integrability (e.g., see [73]) that makes the theory more solvable than one naively expects. We recently 
started exploring the possibility of such structure in A/" = 1 SQCD in [74]. 
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and current two-point function coefficients can be mapped to limitations on the strength 
of gravitational or gauge forces in the presence of light bulk excitations. In the present 
work, we have obtained bounds that scale with the sizes of global symmetry representations, 
which in AdS corresponds to scaling with the size of the bulk gauge group. While many 
of our bounds necessarily apply in a highly quantum regime, we have seen that there 
are at least some bounds (e.g., bounds on operator dimensions in SCFTs) that constrain 
deviations from the large- iV factorization limit, where an AdS description would be weakly 
coupled. It would then be good to find alternate ways of arriving at these bounds in the 
context of AdS, particularly since these constraints are not obvious from the perspective of 
effective field theory [56]. One hopes that thinking more along these lines will lead to a 
deeper understanding of which low-energy theories may admit consistent UV completions, 
particularly in the context of quantum gravity. 



Acknowledgements 

We thank Nima Arkani-Hamed, Diego Hofman, Ken Intriligator, Juan Maldacena, Riccardo 
Rattazzi, Slava Rychkov, and Matt Strassler for helpful comments and conversations. The 
computations in this paper were run on the Odyssey cluster supported by the FAS Science 
Division Research Computing Group at Harvard University. We would like to thank John 
Brunelle in particular for technical support. This work is supported in part by the Harvard 
Center for the Fundamental Laws of Nature, NSF grant PHY-0556111, the Swiss National 
Science Foundation under contract No. 200021-125237, and by the Director, Office of Science, 
Office of High Energy and Nuclear Physics, of the US Department of Energy under Contract 
DE-AC02-05CH11231. 



A Polynomial Approximation Details 

In this appendix we give further details of our implementation of the optimization problem 
discussed in section 2.5 using semidefinite programming. In all of the situations we consider, 
the problem is to find the optimal set of coefficients amnk-, which minimizes the combination 

— amnfcKf^'"*'^'^(0), subject to the constraints 

a„„fc<«''""'(Ao) = 1, (A.l) 
«'mnfcV^^'"*'^''(A) > 0, for all other (non-unit) operators in the spectrum. (A. 2) 

Here V^^''""'^(A) = d'^d'!^Vj^\ denotes derivatives of the fc-th component of the appropriate 

vector £, which may be any of the functions {-Fa,^, -f^A,^, -^a/, -^a/, ^a,^}- The index 
/ denotes possible global symmetry representations. 

As discussed in section 2.5, to apply semidefinite programming we must approximate 
Vf^''""'^(A) as X£(A)P/'"*'^'^(A), where x^(A) is a strictly positive function, and P/''""'^(A) is 
a polynomial in A. Let us begin by discussing derivatives of Fa/ and H/^^i. It is convenient 
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to first rescale each of tliese by a (A, 
sums of terms that factorize: 



-independent function of z and so that they become 
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where ki3{z) = 2;'^/^2-^i(/3/2, /3/2, /3, z). Derivatives of these quantities at (1/2, 1/2) can then 
be straightforwardly evaluated using [52] 

d:[z'-'+^/\E,{P/2,P/2,^,z)]^^^^^ 
T{f3/2 + 2-d) 



a 



2n+{d-l)-/3/2_ 



n] 



5F2(/3/2 + 2 - c/, 13/2, (3/2- 13/2 + 2 - d - n, I3; 1/2) 



2d^ + 8d-2) C^-/ 



T{l3/2 + 2-d 
= 2(5 -2d- n)C^-/ + 2 (/3(/3 - 2) + 2n{n - 3) 
+ 8in-2){n + d-AfCJ^;/ 

= pmkp{i/'i)+QmK{^m- (a.5) 

Here P^{P) and Q^{(3) are polynomials in /3 that can be determined through the above 
recursion relation for Note that taking z ^ 1 — z simply introduces an overall factor 

off-1)'^. 



In Eq. (A.5), we have written derivatives of kp{z) at z = 1/2 in terms of polynomials in 
(3, up to two non-polynomial quantities: /3ki^{l/2) and /c^(l/2). For the purposes of writing 
positivity constraints, we are free to divide by fc^(l/2)//3, which is positive for all (3 that 
occur in unitary theories (/3 > —1). Now, the crucial fact for us is that the remaining 
non-polynomial quantity Kjs = /3/c/3(l/2)//c^(l/2) is meromorphic in (3, and admits a simple 
approximation in terms of rational functions 



Pkp{l/2) 
k'^{l/2) 



M 

4.n 



it 



Nm{I3) 
Dm{I3)' 



(A.6) 



where NM{f3) and Dm{P) are polynomials in /3 of degree M + 1. Here, Vj is the j'th zero of 
[3kp{l/2) and Sj is the j'th zero of /c^(l/2), both of which are close to —2j — 1. Ordinarily 
we would need to account for both the zeros and poles of l3kp{l/2) and ^^(1/2) in the above 
product representation. However, the poles of [3kj3{l/2) and A;^(l/2) coincide at the negative 
odd integers, and so cancel between numerator and denominator.^^ 



^^This recursion relation follows from the hypergeometric differential equation for fc^(z), which itself is a 
consequence of the fact that conformal blocks are eigenfunctions of the quadratic Casimir of the conformal 
group. 

^^Thc factor l/\/2 is lim^_j.oo Kp (with an arbitrary phase for /?), as can be verified using the standard 
integral formula for 2^1 hypergeometric functions. Since this limit exists, Kp is meromorphic on the Riemann 
sphere, not just C. 
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The approximation Eq. (A. 6) becomes arbitrarily good as more zeros are included, and 
moreover converges very quickly. In fact, one can show that 



■l-2j + 0(2-5-5^) J = 0,1,2,..., (A.7) 



so that 




Consequently, it is sufficient to take M ~ a few to achieve an accurate rational approximation 
for that holds uniformly for all physical values /3 > — 1. In practice, we found that 
M = 3 or 4 gives excellent results, which remain effectively unchanged when M is increased. 
Henceforth, we will assume that some appropriate M has been chosen, and write simply 
N{(3) and D{(3) for brevity. 

Combining Eqs. (A. 5) and (A. 6), we can now write 

A;' (1/2) 

where 



u^M ^ N{(3)P^{(3) + HDWm (A.IO) 

is a polynomial in /3, and it can be verified that the pre-factor k'i^{l/2)/ PD{/3) is positive for 
all f3 > —1. Note that the degree of u^{(3) depends on the number of roots M + 1 included 
in the approximation of Eq. (A. 6). 

Derivatives of E^^i^±{z,z) at (1/2, 1/2) can now be written 

9r5^i^A/,±(i/2,i/2) = x.(A)f/ri:±(A), (a.h) 

where 

(A) = fek+.(l/2)fc^„,„,(l/2) 

' ~ {A + £){A-i-2)D{A + £)D{A-£-2) ^ ' ' 

is positive, and 



[/,7i(A) = (l±(-l)-+")K(A + £)«^(A-£-2)-(mOn)] (A.13) 

is a polynomial in A. The inequalities amnkVi {^) > given in (A. 2) are then equivalent 
to a set of polynomial inequalities, which can be rewritten in terms of a semidefinite program 
as described in section 2.5. 

Next let us consider derivatives of the functions J^a,£{z,z) and 'Ha/(-2, ^), appearing 
in superconformal crossing relations. We can again take derivatives using Eq. (A. 9) after 
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rescaling by the same functions of z and z appearing in Eqs. (A. 3) and (A. 4). Applying 
at (1/2, 1/2) to the resulting functions gives 



Tjran 



(A) 



(A 



4(A + ^ + l)D(A + £- 
(A-£-2) D(A-£-2) 



4(A-£-l) D(A- 
(A + £)(A-£-2) 



-/CA+,f/,7i,,,i(A+l) 

^A-£-2 t^rid.ilA + 1) 



D(A + £)Z}(A-£-2) 



16(A + £ + 1)(A - £ - 1) D(A + £ + 2)D(A - £) 



/Ca+£/Ca-£-2 LCT'j_(A 



where 



/3 fc;^+2(i/2) 

/3 + 2 ^^,(1/2) ' 



2) 
(A.14) 

(A.15) 



We can then use the fact that /C^ can be arbitrarily well approximated by a rational function 



/C/3 ~ (12-8V2 



(/3 + l)n.(/3 + 2-..) 



(A.16) 



Again, the approximation improves as more roots are included, and converges after only a 
few terms. Thus, by isolating the polynomial numerator and denominator of the quantity 



/3 



Dm 



-ICf- 



4(/3 + l)D(/3 + 2)'^^- 
we can write the derivatives as a positive function times a polynomial in A: 
XKA) 



(A.17) 



P(A + £)P(A-£-2) 



[ P(A + £)P(A -1-2) f/;7±(A) 



-Ar(A + £)P(A 
^ I?( A + £)A/'( A 
-Ar(A + £)Ar(A 



^-2)^rM,±(A + l) 
^-2)t/rM,±(A + l) 
£-2)f/-^(A + 2)]. 



(A.18) 



Finally, let us note that the results for 2) and HA/iz,!) are identical, but with 

odd-spin terms having the opposite sign. Thus, we see that we can reformulate any of the 
sum rules appearing in SCFTs as a semidefinite program, following the logic described in 
section 2.5. 



B Implementation in SDPA-GMP 

In this appendix we'll give further details of our implementation of the SDP. As we described 
in section 2.5 and appendix A, the general problem (phrased as a SDP) is to minimize 
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—aiVQ^'\0), subject to the constraints 
a.V,y{Ao) = 1, 

a,P/'^(A,(l + x)) = [x]lAl[xU + for < £ < L, 

A^,5/^0 forO<£<L. (B.l) 

For brevity we here we use the index i = 1, . . . , k{k + l)/2 x dim(y^^) to run over all of 
the z and z derivatives under consideration, as well as the components of the vector V^^ ^. / 
runs over possible global symmetry representations, and and Bj, are positive semidefinite 
matrices. We recall that [x]d is the vector with entries . . . and if the polynomial 
P/'* has degree 27^ + 1 — (with = 0, 1), then di = 27^ and d'^ = 27^ — 2e£. 

The middle constraint is an equality between polynomials in x, so in practice we will 
implement it by matching each polynomial coefficient: 



= coeffs,. 



-a,P/'X A,(l + x)) + TTiXaXx)Al) + xTt{XMB',)] . (B.2 



In this expression we have also defined the matrix Xd{x) = J. Since many SDP solvers 

only allow positive variables, in practice it will additionally be convenient to introduce a 
'slack variable' s, where without loss of generality we can replace — >■ — s in the above 
expressions and require a^, s > 0. 

We solve the above semidefinite program using SDPA-GMP 7.1.2 [82], which utihzes the 
GNU Multiple Precision Arithmetic Library (GMP). We use Mathematica 7.0 to compute 
the vectors V/'^ and polynomials P/'*, performing all computations using 100 digits of 
precision. When using the approximations of Eqs. (A. 6) and (A. 16) we keep four roots, 
leading to approximations that differ from the exact functions by ~ 10~^ — 10~^*^, depending 
on the value of (3. In our computations we have found it sufficient to take L = 20; in 
addition we add constraints for i = 1000, 1001 in order to effectively include the asymptotic 
constraints at large i. After setting up the problem in Mathematica, we write the SDP to a 
file using the SDPA sparse data format. 

When running SDPA-GMP, we use the parameters: 



SDPA-GMP Parameter 


Value 


maxlteration 


1000 


epsilonStar 


10-10 


lambdaStar 


1020 


omegaStar 


1020 


lowerBound 


_1020 


upperBound 


1020 


betaStar 


0.1 


betaBar 


0.3 


gammaStar 


0.9 


epsilonDash 


10-10 


precision 


200 
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To make our plots, we run data points in parallel using the Odyssey computing cluster at 
Harvard University. In the majority of our plots we use a horizontal spacing of 5d = 10"^, 
supplemented by a higher resolution scan with 6d = 10~^ for d < 1.01 {6d — )■ ^Aq in 
figure 10). To compute dimension bounds, we vary Aq using a binary search, terminating 
at a vertical resolution of 10~^. In all cases that we have checked, increasing L or including 
more roots in the polynomial approximation leads to a completely negligible (< 10~^) change 
in the computed bound. 
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